AD-A274  586 

■llllllll 


WL-TR-93-4098 

OPTICAL  CHARACTERIZATION  AND  GUIDED 
WAVE  LOSS  MEASUREMENTS  FOR  END-CAPPED 
6F-POLYBENZOXAZOLE 


ANGELA  L.  MCPHERSON 
JEFFERY  W.  BAUR 


SEPTEMBER  1993 

FINAL  REPORT  FOR  11/01/90-06/01/93 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


f 


MATERIALS  DIRECTORATE 
WRIGHT  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AFB  OH  45433-7734 


94  1  7  315 


NOTICE 


WHEN  GOVERNMENT  DRAWINGS.  SPECIFICATIONS,  OR  OTHER  DATA  ARE  USED  FOR 
ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFINITELY  GOVERNMENT- 
RELATED  PROCUREMENT,  THE  UNITED  STATES  GOVERNMENT  INCURS  NO  RESPONS¬ 
IBILITY  OR  ANY  OBLIGATION  WHATSOEVER.  THE  FACT  THAT  THE  GOVERNMENT  MAY 
HAVE  FORMULATED  OR  IN  ANY  WAY  SUPPLIED  THE  SAID  DRAWINGS,  SPECIHCATIONS, 
OR  OTHER  DATA,  IS  NOT  TO  BE  REGARDED  BY  IMPLICATION,  OR  OTHERWISE  IN  ANY 
MANNER  CONSTRUED,  AS  LICENSING  THE  HOLDER,  OR  ANY  OTHER  PERSON  OR 
CORPORATION;  OR  AS  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE,  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED 
THERETO. 


THIS  REPORT  IS  RELEASABLE  TO  THE  NATIONAL  TECHNICAL  INFORMATION  SERVICE 
(NTIS).  AT  NTIS,  FT  WILL  BE  AVAILABLE  TO  THE  GENERAL  PUBLIC,  INCLUDING  FOREIGN 
NATIONS. 


THIS  TECHNICAL  REPORT  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR  PUBLICATION. 


Ang^a  L.  McPherson,  Research  Physicist 
Electronic  &  Optical  Mat'ls  Branch 
Eleccrotnagnedc  Mat'ls  &  Surv.  Div. 


ROBERT  L.  DENISON,  Chief 
Electronic  &  Optical  Mat’ls  Branch 
Electromagnetic  Mat’ls  &  Surv.  Div. 


WILLIAM  R.  WOODY,  ^ief 
Electromagnetic  Mat'ls  &  Surv.  Div. 
Materials  Directorate 


IF  YOUR  ADDRESS  HAS  CHANGED,  IF  YOU  WISH  TO  BE  REMOVED  FROM  OUR  MAILING 
LIST,  OR  IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY  YOUR  ORGANIZATION  PLEASE 
NOTIFY  WL/MLPO.  BLDG  651,  3005  P  ST  STE  6,  WPAFB  OH  45433-7707,  TO  HELP  MAINTAIN 
A  CURRENT  MAILING  LIST. 

COPIES  OF  THIS  REPORT  SHOULD  NOT  BE  RETURNED  UNLESS  RETURN  IS  REQUIRED 
BY  SECURITY  CONSIDERATIONS,  CONTRACTUAL  OBLIGATIONS,  OR  NOTICE  ON  A 
SPECinC  DOCUMENT. 


REPORT  DOCUMENTATION  PAGE 

form  Approved 

0MB  No.  0704-0188 

Public  rtpeiting  burdpn  for  this  colitction  of  informotion  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sOMroes. 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  colicction  of  information.  Send  comments  regarding  th^  burden  estimate  or  any  otiier  aspect  of  this 
colliKtion  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Direaorate  for  information  Operations  and  Reports.  121S  Jefferson 
O^is  Highway.  Suite  1204.  Arlington.  VA  22202*4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Proiect  (0704-0 18B).  Washington.  OC  20S03. 

1.  A€ENCY  USE  ONLY  (Lejve  blank)  2.  REPORT  OATE  3.  REPORT  TYPE  AND  OATES  COVERED 

SEP  1993  FINAL  11/01/90—06/01/93 

A.  TITLE  AND  SUBTITLE  OPTICAL  CHARACTERlZATlOU  AND  C3UIDED 
WAVE  LOSS  MEASUREMENTS  FOR  END-CAPPED  6F- 
POLYBENZOXAZOLE 

S.  FUNDING  NUMBERS 

PE  61102 

PR  2305 

TA  FO 

WU  FJ 

«.  AUTHOR(S) 

ANGELA  L.  MCPHERSON 

JEFFERY  W.  BAUR 

7.  PERFORMING  ORGANIZATION  NAME(S}  ANO  ADORESS(ES} 

MATERIALS  DIRECTORATE 

WRIGHT  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT  PATTERSON  APB  OH  45433-7734 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

WL-TR-93-4098 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 

MATERIALS  DIRECTORATE 

WRIGHT  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRII3HT  PATTERSON  APB  OH  45433-7734 

10.  SPONSORING /MONirORMG 

AGENCY  REPORT  NUMBER 

WL-TR-93-4098 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS 
UNLIMITED. 

12b.  distribution  CODE 

13.  ABSTRACT  (Maximum  200  words) 

A.  Mathematica  software  package  ±0.  modeJ  . the  i>ehavrior  of.  light  inside  an  optical 
waveguide  is  discussed.  Examples  of  package  features  are  presented  and  discussed. 
Possible  new  waveguide  materials  are  investigated  with  this  package  as  a  baseline 
for  material  evaluation.  The  use  of  the  polymer  6F-PB0-1E  as  a  possible 
candidate  for  high  temperature  optical  waveguide  material  based  in  its  high  Tg  of 
300C  and  its  optical  clarity  is  discussed.  Theoretical  studies  are  performed 
with  the  Mathematica  package.  Slab  waveguide  samples  are  fabricated  by  spin¬ 
coating  on  oxidizOS  silicon  wafers  and  characterizOd  through  measurements  of  the 
index  of  refraction  and  waveguide  loss.  Results  show  that  the  material  exhibits 
its  lowest  loss  of  4.5  dB/cm  at  a  wavelength  of  700nm.  We  conclude,  based  on 
experimental  results,  that  this  material,  6F-PB0-1E,  while  it  may  be  useful  for 
some  optical  applications,  does  not  meet  the  Air  Force  specification  of  1  dB/cm 
or  less  performance  guideline  for  optical  waveguide  materials  and  therefore  is 
not  appropriate  for  Air  Force  applications  in  its  present  form. 

14.  SUBJEa  TERMS 

optical  waveguide  loss,  planar  waveguide,  absorption,  index 
of  refraction 

IS.  NUMBER  OF  PAGES 

96 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRAa 

_ Ub _ 

MSN  7S40-01-280-S500  Standard  Form  298  (Rev  2-89) 


Pretcnbed  by  ANSI  Std  239-)8 


2W-102 


TABLE  OF  CONTENTS 


Section  Page 

1.0  INTRODUCTION  1 

2.0  THEOREnCAL  BACKGROUND  2 

2.1  WAVEGUIDE  PARAMETERS  2 

2.2  EIGENVALUE  EQUATION  3 

2.3  LIMITS  ON  n«  6 

3.0  MODELING  9 

3.1  PACKAGE  FEATURES  11 

3.1.1  PROCEDURE  "diecklnput’*  11. 

3.1.2  PROCEDURE  ■•PlotThis"  12 

3.1.3  PROCEDURE  "FindRoots"  12 

3.1.4  PROCEDURE  "PrppagationConstant"  13 

3.1.5  ATTENUATION  CALCULATION  PROCEDURES  16 

3. 1.5.1  PROCEDURE  "N1  Attenuation"  16 

3. 1.5.2  PROCEDURE  "N2Attenuation"  18 

3. 1.5.3  PROCEDURE  "N3Attenuation"  19 

3. 1.5.4  PROCEDURE  "N4Attenuation"  20 

3. 1.5.5  PROCEDURE  "GuideAttenuation"  21 

3. 1.5.6  PROCEDURE  "OaddingAttenuation"  22 

3. 1.5.7  PROCEDURE  "WavelengthAttenuation"  23 

3.1.6  PROCEDURE  "CutoffThickness"  24 

3.1.7  P  :EDURE  "CutoffWavelength"  25 

3.2  MODELING  EXAMPLE  26 

iii 


TABLE  OF  CONTENTS  {Continued) 

Section  Page 

4.0  WAVEGUIDE  MATERIAL  30 

4.1  WAVEGUIDE  FABRICATION  31 

4.1.1  THE  SUBSTRATE  31 

4.1.2  THE  GUIDING  MATERIAL  35 

4.2  SAMPLE  CHARACIERIZATION  36 

4.2.1  GUIDING  LAYER  THICKNESS  36 

4.2.2  INDEX  OF  REFRACTION  37 

4.2.3  ABSORPTION  40 

4.3  THEORETICAL  BEHAVIOR  43 

5.0  EXPERIMENTAL  MEASUREMENTS  OF  WAVEGUIDE 

ATTENUATION  44 

5.1  INTRODUCnON  44 

5.2  LOSS  MEASUREMENT  TECHNIQUES  45 

5.2. 1  THE  THREE-PRISM  METHOD  45 

5.2.2  THE  TWaPRISM  METHOD  47 

5.2.3  THE  OUT-OF-PLANE  SCATTERING  METHOD  47 

5.3  EXPERIMENTAL  PARAMETERS  49 

5.4  RAW  WAVEGUIDE  DATA  50 

5.5  THREE  METHODS  COMPARISON  53 

5.6  WAVEGUIDE  RESULTS  54 

6.0  CONCLUSION  57 

7.0  REFERENCES  58 

APPENDIX  A:  BASE  CLEAN  PROCEDURE  60 


IV 


Section 


TABLE  OF  CONTENTS  (Concluded) 


Page 


APPENDIX  B:  RAW  WAVEGUIDE  DATA  TABLES  61 

APPENDIX  C:  PlanarWaveguidcs.m  68 

DOCUMENTATION  68 

SOURCE  CODE  69 


DTIG  OfTTALITY  INSPECT  9 


Ml  0^'. 

:  ■  □ 
ih  i  □ 


1  .  CC'’33 


LIST  OF  FIGURES 

Hgure  Page 

1  Diagram  of  the  Four  Layer  Planar  Waveguide  2 

2  Exan^le  of  '‘PtopagationQmstant"  Procedure  IS 

3  Example  of  "NlAttenutdcm"  Procedure  17 

4  Example  of  "NlAttenuatitm''  Procedure  18 

5  Example  of  ”N3Attenuation''  Procedure  19. 

6  Example  of  "N4Attenuation"  Procedure  20 

7  Exan^le  of  "GuideAttenuadon"  Procedure  21 

8  Example  of  "QaddingAttaiuadon"  Procedure  22 

9  Example  of  "WavelengthAttenuadon"  Procedure  23 

10  Example  of  "CutofifThickness"  Procedure  24 

1 1  Example  of  "CutoffWavclength"  Procedure  25 

12  Theoredcal  Calculation  of  Cutoff  Wavelength  for  0.4^m 

PMMA  on  0.65^m  SiC^  on  Si  27 

13  Theoretical  Calculation  of  the  Propagation  Constant  for 

0.4pm  PMMA  on  0.65pm  SiC)2  on  Si  at  A.  =  514.5nm  28 

14  Theoretical  Calculation  of  the  Propagation  Constant  for  0.5pm 

PMMA  on  2.0pm  Si02  on  Si  at  A.  =  5 14.5nm  29 

15  6F-PBO  and  Variations  31 

16  Oxide  Thickness  Profile  taken  on  Profilometer  for  3-Inch 

Oxidized  Silicon  Wafer  33 

17  Refiactive  Index  Profile  of  Oxide  Layer  taken  on  Ellipsometer 

for  3-Inch  Oxidized  Silicon  Wafer  33 


VI 


LIST  OF  FIGURES  (Continued) 

Figure  Page 

18  Oxide  Thickness  Profile  taken  on  Profilonoeter  for  2-Inch 

Oxidized  Silicon  Wafer  34 

19  Refractive  Intkx  Profile  of  Oxide  Layer  taken  on  Ellipstxneter 

fcx*  2-Inch  Oxidized  Silicon  Wafer  34 

20  Polymer  Thickness  Profile  taken  on  Profilometer  for  2-Inch 

Wafer  37 

21  Refractive  Index  Profile  of  Polymer  Layer  taken  on  Ellipstxneter 

for  2-Inch  Wafer  39 

22  Transmission  Data  taken  on  Free-Standing  Him  of  6F-PBO-  IE 

with  Spectrophotometer  41 

23  Absorption  Coefficient  for  6F-PBO'  IE  as  Calculated  frran  the 

Transmission  Data  42 

24  Theoretical  Plot  of  Loss  versus  Wavelength  Computed  with 

Mathematica  43 

25  Experimental  Setup  for  the  Three-Prism  Method  to  Measure 

Attenuation  45 

26  Experimental  Setup  for  the  Out-of-Plane  Scattering  Method, 

Showing  Detector  Positioning  and  Masking  48. 

27  Raw  Data  of  l.O^m  6F-PBO-1E  on  1.4pm  Si02  on  Si  at  488nm  50 


28  Raw  Data  of  1 .0pm  6F-PBO-1E  on  1 .4pm  Si02  on  Si  at  5 14.5nm  5 1 

29  Raw  Data  ofl.Opm6F-PBO-lEon  1.4pm  Si020nSi  at  578.15nm  51 

30  Raw  Data  of  1.0pm  6F-PBO-1E  on  1.4pm  Si02  on  Si  at  800nm  52 

31  Raw  Data  of  1.0pm  6F-PBO-1E  on  1.4pm  Si02  on  Si  at  900nm  52 


Vll 


LIST  OF  FIGURES  {Concluded) 

Figure  Page 

32  Raw  Data  of  l.O^m  6F-PBO-1E  on  1.4)iin  Si02  on  Si  at  975nm  53 

33  Con^arison  of  Three  Mediods  of  Loss  Measurenoent  on  1  .O^m 

6F-PBO-IE  on  1.4pin  Si02  on  3"  Si  Wafer  at  900nin  54 

34  Summary  of  Loss  Data  for  6F-PBO-  IE  at  0.45  -  0.975^m 

Including  Out-of-Plane,  Three-Prism  and  Two-Prism  Methods  55 


Vlll 


List  of  Tables 


Table  Page 

1  Compaiison  of  Three  Methods  for  Finding  Index  of  Refraction 

and  Thickness  of  6F-PBO-  IE  39 

2  Summary  of  Average  Loss  Data  for  6F-PBO-  IE  56 


ACKNOWLEDGEMENTS 


The  authors  wish  to  acknowledge  S  J.  Bai  of  the  Materials  Directorate  for  his  help 
in  verification  of  the  index  of  refraction  measurements  and  Dr.  Uma  Ramabadran  of  the 
University  of  Dayton  Research  Institute  for  her  support  and  consultation.  Mr.  Ronald  O. 
Neaffer  is  acknowledged  for  his  assistance  in  the  ellipsometric  work  and  Mr.  E.  Thomas 
Kensky  for  his  assistance  with  the  waveguide  stage  and  progranoming  efforts.  We  also 
wish  to  acknowledge  Daychem  Inc,  Centerville  OH,  for  supplying  the  polymer  material  on 
which  this  study  is  based.  Finally,  we  wish  to  thank  Dr.  F.  Ken  Hopkins  and  Dr.  David 
Zelmon  for  their  help  and  suggestions  on  the  modeling  section  of  this  report. 


1.0  INTRODUCTION 


Optical  waveguides  are  structures  that  carry  optical  signals  by  confining  light  in  one 
or  two  dimensions  and  have  been  studied  for  some  time*.  These  structures  have  been 
studied  extensively  during  the  past  two  decades  because  they  are  the  central  component  of 
most  integrated  optical  devices  including  couplers,  nnodulators,  and  interferometers.  Given 
the  wide  variety  of  materials  and  processes  that  may  be  used  in  waveguide  fabrication,  it  is 
imperative  that  designers  be  able  to  predict  the  performance  of  this  critical  component.  We 
have  chosen  to  develop  calculations  for  the  four  layer  planar  waveguide  structure  since  it  is 
generally  adequate  to  model  the  behavior  of  light  energy  in  various  device-specific 
waveguide  configurations  such  as  charmel  or  rib-type. 

We  have  developed  a  program  that  may  be  used  to  calculate  field  distribution, 
propagation  loss,  and  critical  values  for  the  four-layer  planar  optical  waveguide.  This 
problem  has  previously  been  solved  using  a  variety  of  methods^.  We  discuss  a  simpler 
approach  to  solving  the  complex  eigenvalue  equation  using  the  software  package 
Mathematica. 
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2.0  THEORETICAL  BACKGROUND 


2.1  THE  WAVEGUIDE  PARAMETERS 

A  planar  waveguide  is  a  multilayer  structure  such  as  that  shown  in  Hgure  1 .  At 
an  interface  between  two  media  of  different  indices  of  refraction,  electromagnetic  (EM) 
waves  are  split  into  a  reflected  and  refracted  beam.  The  reflected  beam  obeys  the  law  of 
reflection,  and  the  refracted  beam  obeys  Snell's  Law : 

Ilisin  01  =  n2sin  02  (1) 

where  ni,  n2  are  the  indices  of  refraction  for  regions  1  and  2  and  0i,  02  are  the  angles  of 
incidence  and  refraction,  respectively,  and  are  defined  with  respect  to  the  normal  to  the 
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interface  between  regions  1  and  2.  In  the  case  ni<n2  there  is  always  a  real  angle.  02,  that 
satisfies  this  equation  and,  therefore,  there  is  transmittance  from  region  1  to  region  2. 
However,  if  ni>n2  and  02>0i,  then  there  is  a  maximum  angle  of  incidence  for  which  there 
is  a  real  solution  for  02.  This  occurs  when  02  reaches  its  maximum  value  of  90®.  The 
corresponding  value  of  the  incident  beam  angle  is  known  as  the  critical  angle  and  is  given 
by; 

^in  0c  =  1  (2) 

Hi 

Total  internal  reflection  occurs  at  angles  of  incidence  larger  than  the  critical  angle  and, 
therefore,  there  is  no  refracted  beam  (02=90®).  If  a  material  with  refractive  index  n2  is 
placed  between  two  other  materials  widi  indices  ni  and  n^  so  that  ni4i3<n2,  the  light  can  be 
confined  to  the  center  layer  by  total  internal  reflection.  This  is  the  situation  in  the  planar 
waveguide. 

The  structure  under  consideration  (see  Figure  1)  consists  of  four  layers,  each 
having  a  different  index  of  refraction.  Light  is  coupled  into  the  waveguide  from  layer  1, 
usually  air.  Layer  2,  called  the  guiding  layer,  is  the  high-index  layer  where  the  energy  is  to 
be  confined.  Layer  3  is  the  buffer  layer  and  is  used  to  prevent  any  possible  coupling  of  the 
guided  wave  from  layer  2  to  layer  4.  The  entire  arrangement  is  fabricated  on  a  substrate, 
layer  4,  which  can  have  a  very  high  index. 

2.2  THE  EIGENVALUE  EQUATION 

The  behavior  of  light  in  a  structure  such  as  Rguie  1  is  calculated  by  solving  the  wave 
equation  in  each  of  the  four  layers  and  applying  appropriate  boundary  conditions.  In  each 
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layer  of  the  gtiide,  the  behavior  of  the  EM  wave  E(x,y,z,t)  is  described  by  the  wave 
equation^ ; 


V^f)  +  kin2i(r)  =  0  (3) 

where  ko=27C/X.  X  is  the  wavelength  of  the  propagating  light,  and  n  is  the  refractive  index 
of  the  material.  We  assume  the  waves  are  time  harmonic,  the  dielectric  is  lossless  and  there 
are  no  free  charges. 

The  waveguide  problem  may  be  solved  by  assuming  that  the  optical  fields 
propagate  as  inhomt^eneous  plane  waves  of  the  fcmn: 

S(r.t)=El(x.y)e+i<‘“-P*)  W 

Assuming  no  y-ditection  variation,  substimting  this  expression  into  the  wave  equation 
yields: 

1^1  eJx)  +  [k2  n?  -  eJx)  =  0  (5) 

where  ni  is  the  index  of  the  i-th  layer 

We  define: 


to  be  the  propagation  parameters  of  layars  1, 2, 3,  and  4,  respectively. 

(Note:  P  =  ko  ne  where  n*  is  the  unknown  effective  refractive  index  of  the  waveguide.) 

The  wave  equation  will  have  two  sets  of  solutions  coiresponding  to  the  two  orthogonal 

polarizations:  transverse  electric  (TE)  and  transverse  magnetic  (TM).  Calculations  for  both 
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TE  and  TM  modal  solutions  proceed  in  the  same  manner,  so  our  discussion  will  be 
confined  to  die  TE  modes  of  prcqiagation  only. 

For  the  four-layer  guide,  we  assume  the  relation  n2>n3>ni.  Remembering  this  and 
going  back  to  equation  (S),  we  can  write  down  the  solutions  to  the  equation  in  each  layer  of 
the  guide.  The  desired  condition  consists  of  en^gy  confinement/propagation  in  the  guiding 
layer,  evanescent  fields  in  the  first  and  thinl  layers  and  possible  propagation  in  the  fourth 
layer.  We  can  solve  the  ordinary  differential  equations  (ODE)  and,  realizing  that  all  fields 
are  finite  and  no  fields  originate  in  the  substrate,  the  solutions  are: 

^(r,t)  =  Ai{x,y)e<*-*) 

Ejr,t)  =  A^x,y)  sin  (hx)  +  A2'(x,y)  cos  (hx)  (7) 

E3(r,t)  =  A3{x,y)  eP<*+*)  +  A3'(x,y)  c-p<*'“> 

E4(r,t)  =  A4x,y)  c'‘  *(*+*<*) 

Now  we  must  impose  boundary  conditions.  At  each  interface,  in  the  absence  of  any 
sources,  both  the  tangential  electric  (E)  and  tangential  magnetic  (H)  fields  must  be 
continuous.  For  source-fiee  media,  the  E  and  H  fields  are  related  via  Maxwell's  equatitxis: 

VxE  =  -icojiH 
Vxii  =  itoeE 


The  cotiesponding  H  fields  are: 


Haff.t)  =  [Az  h  cos  (hx)  -  A2  h  sin  (hx)]  (9) 

iSf  .t)  =  ^  [  As  p  eK”*)  -  A,  p  e-P(»«)  ] 
S(M=7i7[A4Se-is9‘«*'»] 
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Upon  application  of  the  boundary  conditions,  the  following  set  of  equations  result  (where 
a,  b  are  defined  in  Figure  1). 

Ai  =  Ai  sin  (ha)  +  A2  cos  (ha) 

A2  sin  (-ha)  +  A2  cos  (-ha)  ==  A3  +  A3 

K  •  w  (10) 

A3e-i*+ A3eP‘>=  A4 

-Ai  q  =  A2  h  cos  (ha)  +  A2  h  sin  (ha) 

A2  h  cos  (-ha)  -  A2  h  sin  (-ha)  =  A3P- A3P 
i  [A3  p  e'P**  -  A3  p  eP‘’]=  s  A4 


This  system  of  equations  can  be  simplified  to  one  equation  with  one  unknown,  this  P 
may  be  complex,  representing  energy  lost  to  the  system  via  coupling  into  the  substrate. 
The  eigenvalue  equation  for  ^  is: 

(11) 

(p  +  is)  [(h^  -  pq)  Tan(2ah)  -  h(p  +  q)]  +  e-2pt  (p  -  is)  [(h^  +  pq)  Tan(2ah)  +  h(p  -  q)l  =  0 
where  t  =  (b-a), 

n4  =  n4real  ■*"  i  tHimag, 

P  =  Preal  1  Pimag 

2.3.  LIMITS  ON  n® 

The  limits  on  the  complex  propagation  constant,  P,  must  be  examined  before  any 
modeling  is  presented.  These  limits  form  a  baseline  and  are  used  extensively  in  arriving  at 
the  soluticHi  to  the  eigenvalue  equation. 

It  has  been  shown  that  to  enable  guided  mode  solutions,  the  index  of  redaction  of 
the  guiding  layer  must  be  higher  than  that  of  the  surrounding  layers.  The  assumption  is 
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also  made  that  n2>n3>ni,  which  is  always  true  if  the  first  layer  is  air.  Next,  we  recall  the 
wave  equation  in  the  general  form  inside  each  layer  as  shown  in  equation  (S).  It  has  the 


clause: 

[]^n?-P^]  (12) 

which  will  define  the  behavior  of  the  equation.  In  the  first  case,  take  P  >  kon2.  Simplifying 
die  equatitxi  above  results  in  the  following  generalized  ODE: 


i—]>o 

El  3x2  I 


This  can  be  rewritten  as: 


(13) 


32 

.3x2  I 


=  a2E 


(14) 


where  a  is  some  constant  The  general  solution  to  this  ODE  is  of  the  form: 


E  =  Ae«*  +  Be-«x  (15) 

This  solution  then  corresponds  to  all  layers  in  the  structure  based  on  the  assumptitm  that 
n2>n3>ni.  Reviewing  the  defmititm  of  propagation  as  outlined  in  Section  2.1,  diis  solution 
does  not  correspond  to  a  guided  mode. 

The  second  case  investigated  is  kon3  <  P  <  kon2.  In  this  case,  the  ODE  becomes: 


3x2. 


=  -  a^E 


(16) 


The  two  regions,  P  <  kon2  and  koni,kon3  <  P  must  now  be  examined.  P  now  has  two 
solutions  depending  on  which  region  IS  examined.  Inlaym'2, 

E  =  Asin(ax)  +  B  cos(ax)  (17) 

while  in  layers  1  and  3  the  same  solution  as  in  equation  (14)  exists.  This  solution 

corresponds  to  a  guided  mode,  i.e.,  there  is  propagation  in  tbr  guiding  layer  (layer  2),  and 

evanescent  fields  in  the  suiroundircr  layers  (layers  1  and  3). 
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In  the  case  where  koni  <  P  <  kons  there  is  exponential  behavior  in  layer  1  and 
sinusoidal  behavior  in  regions  2  and  3.  Also,  if  0  <  ^  <  koni,  there  is  sinusoidal  behavior 
in  all  regions.  Neither  of  these  corresponds  to  a  guided  mode  solution.  Therefore,  the 
limits  on  p  for  guided  mode  solutions  is  konj  <  ^  <  kon2- 


8 


3.0.  MODELING 


The  eigenvalue  equation  may  be  solved  using  the  mathematics  software  package 
Mathematica  to  calculate  the  field  distributions,  propagation  loss,  and  critical  values  for  the 
four-layer  planar  waveguide.  The  basic  algorithm  calculates  all  possible  roots  of  the 
eigenvalue  equation,  defines  all  confined  modes  and  plots  the  waveform  for  each  mode. 
The  program  requires  only  specification  of  the  refnctive  indices  of  the  four  layers,  the 
thickness  of  the  guiding  and  buffer  layers  and  the  wavelength  of  the  input  light 

This  algorithm  can  also  determine  other  waveguide  pn^>erties  including  waveguide 
attenuation  (energy  loss  into  the  fourth  layer)  and  the  cutoff  values  (wavelength  and 
guiding  layer  thickness  below  which  no  modes  are  confined)  which  are  important 
ctxisiderations  in  fatnication. 

Mathematica'%  built-in  root-rinding  procedure,  FindRoot,  is  used  to  solve  the 
equation.  The  solution  must  lie  between  the  values  k  n3  and  k  n2  (k  is  ItUX,  X  being  the 
wavelengdi  of  the  incident  light)  for  the  equation  to  hold  true.  We  can  use  this  stimulation 
to  define  an  initial  value  for  FindRoot,  however,  the  root  converges  very  slowly.  A 
quieter  method  is  to  solve  the  problem  by  considering  an  infinitely  duck  buffer  layer.  This 
is  equivalent  to  considering  a  three-layer  asymmetric  waveguide^  The  four-layer 
eigenvalue  equation  reduces  m  (s  corresptxids  to  the  fourth  layer  and  is  set  to  zero): 

Tan(2ah)=^^  (18) 

Since  this  equation  has  roots  that  converge  relatively  quickly,  the  value  of  P  which  is  the 
solution  to  this  eigenvalue  equation  is  used  as  the  starting  value  for  the  FindRoot  routine 
in  calculating  the  eigenvalues  for  the  four-layer  equation. 
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For  a  given  set  of  parameters,  the  four-layer  eigenvalue  equation  may  have  multiple 
roots  corresponding  to  the  different  modes  supported  by  the  waveguide.  It  is  easy  to  see 
bxnn  the  three-layer  eigenvalue  equadtm  that  a  simple  Do  loop  changing  the  initial  values 
for  FIndRoot  will  not  solve  for  all  roots.  The  equation  blows  up  at  certain  initial  values. 
For  this  reason,  the  combination  of  Mathematica's  FindRoot,  Check,  and  Continue 
functions  are  indispensable.  Check  is  a  boolean  variable  which  scans  the  output  stream 
for  error  messages.  If  errors  are  not  encountered.  Check  remms  the  value  of  its  first 
argumenL  The  value  of  the  second  argument,  in  this  case  Continue,  is  returned  if  errors 
are  found.  Combining  these  functions  (alcxig  with  checks  for  duplicate  roots  and  assuming 
FourLayerEquation  is  the  eigenvalue  equation)  gives  the  root-finding  portion  of  the  code: 

Check[rule2  =  FindRoot[FourLayerEquation  ss  0, 
{B,B/«i‘ulel},AccuracyGoal  >>  12],Continue[]]; 

Check  and  Continue  allow  for  the  evaluation  at  these  points  without  crashing  the  code. 

The  starting  point  'BArulel'  is  the  solution  to  the  much  simpler  three-layer  problem. 

The  source  code  shown  above  will  completely  solve  the  eigenvalue  equation  for  all 

confined  modes  of  the  user-specified  waveguide  structure.  Once  the  propagation  constant, 

P,  is  known,  it  is  used  to  calculate  the  amplitudes  Ai  in  each  layer  of  the  waveguide. 

The  basic  algorithm  scans  the  range  of  possible  roots  of  the  eigenvalue  equation.  It  then 
produces  all  roots,  thereby,  defiiung  all  confined  modes,  and  plots  the  waveform  for  each 
mode.  The  program  requires  only  specification  of  the  refractive  indices  of  the  four  layers, 
the  thickness  of  the  guiding  and  buffer  layers,  and  the  wavelength  of  the  input  light 

The  relative  simplicity  of  this  algorithm  allows  for  the  determination  of  other 
waveguide  properties.  Also  offered  in  this  package  are  routines  to  calculate  the  waveguide 
attenuation  (energy  loss  into  the  fourth  layer  and  associated  with  the  specified  structure) 
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and  the  cutoff  values  (wavelength  and  guiding  layer  thickness  below  which  no  modes  are 
ccMifined)  which  are  important  considerations  in  the  fabrication  of  waveguiding  structures. 

3.1.  PACKAGE  FEATURES 

The  following  subsections  describe  and  give  examples  of  the  procedures  featured  in 
the  Mathematica  package  PlanarWaveguides.m.  For  further  details  on  the  code,  see 
Section  10.0. 

3.1.1.  PROCEDURE  "Checklnput" 

Note:  This  is  an  internal  procedure. 

Several  checks  are  made  on  the  input  values  to  see  if  there  are  any  errors  in  the 
input  which  would  preclude  the  existence  of  any  guided  modes.  These  checks  include 
cutoff  thickness,  cutoff  wavelength,  ni<Re[n2]>n3,  and  other  erroneous  input  values.  If 
any  of  these  errors  are  encountered,  execution  of  the  code  will  cease.  A  default  of  one 
iteration  through  the  root-finding  procedure  is  assumed  for  the  case  where  the  increment 
entered  is  zero. 

The  cutoff  thickness  and  cutoff  wavelength  are  calculated  from  the  three-layer 
eigenvalue  equation  as  the  value  of  the  effective  index  of  refraction  (proportional  to  the 
propagation  constant)  approaches  the  index  of  the  third  (cladding)  layer.  If  the  waveguide 
thickness  is  less  than  the  cutoff  thickness  or  the  cutoff  wavelength  is  shorter  than  the  input 
wavelength,  an  appropriate  message  is  printed  and  program  execution  stops. 


3.1.2.  PROCEDURE  "PlotThis” 

Note:  This  is  an  internal  procedure. 

As  the  number  of  guided  modes  may  differ  as  the  parameters  are  varied,  the  values 
of  attenuation  are  separated  by  mode  number.  "PlotThis”  accepts  the  list  of  roots  found 
and  marks  the  associated  position  in  "modeFlags"  as  to  the  noode  number.  "ModeFlags"  is 
then  used  to  plot  the  attenuation  versus  the  mode  number.  The  Re[newList[[i,l]]]  is  used 
for  the  cases  that  involve  the  complex  na  parameter. 

A  variation  of  tu  (or  n^)  may  result  in  values  such  as  3.8S  -  0.01 1  and  3.85  -  0.02 1. 
The  attenuation  values  are  considered  different  for  the  number  of  modes,  but  the  values  will 
be  ploned  above/below  each  other  as  only  the  Re[n4]  (Re[n2])  is  used  in  plotting. 

The  parameters  for  this  routine  are  tiie  list  of  waveguide  parameter  values  with  their 
associated  attenuation  values,  the  "plotLabel,"  and  "biggestListCount"  which  holds  the 
maximum  number  of  guided  modes  fOT  any  given  value  of  the  varied  parameter. 

3.1.3.  PROCEDURE  "FindRoots” 

Note:  This  is  an  internal  procedure. 

"HndRoots"  is  the  core  of  the  code.  The  roots  of  the  eigenvalue  equation  are 
acmally  found  in  this  procedure. 

Since  the  interest  is  only  in  the  guided  mode  solutions  (roots  of  the  eigenvalue 
equation  between  the  index  of  layer  two  and  three  times  a  constant  k),  the  testing  for  roots 
scans  the  entire  range  of  possible  roots,  with  a  step-size  (increment2)  dependent  on  the  size 
of  the  range. 


The  parameters  are  the  increment  size  for  the  test  range  (increment2),  the  present 
value  of  the  varied  parameter  (variable),  "threeLayerEquation,''  "fourLayerEquadon,"  k, 
n2,  ns,  and  flag  which  is  simply  used  to  test  if  the  guiding  layer  thickness  is  the  parameter 
being  varied.  This  must  be  tested  for  because  the  equadons  and  theory  depend  on  the  value 
"halfGuide"  which  is  one-half  the  guiding  layer  thickness.  To  avoid  confusing  the  user, 
the  value  of  the  total  guiding  layer  thickness  is  printed  out  instead  of  "halfGuide." 

To  obtain  starting  values  for  the  four-layer  algorithm,  the  soludons  for  the  three-layer 
waveguide  are  calculated.  The  simplicity  of  the  three-layer  eigenvalue  equadon  allows  less 
stringent  requirements  on  the  starting  values  for  its  soludons.  These  soludons  can  then  in 
turn  be  used  as  initial  condidons  for  the  four-layer  equadon. 

The  "AccuracyGoal"  and  "Maxiteradons"  were  chosen  to  allow  quick  tum-around-dme, 
yet  sdll  achieve  accurate  soludons.  To  aid  in  shortening  the  calculadon  dme,  if  the  root 
found  is  greater  than  the  test  value,  the  test  value  will  be  set  to  the  new  root  "FindRoots" 
returns  a  list  of  P  values  which  is  added  (Union)  to  the  list  of  all  soludons  found  for  the 
varied  parameter. 

"Bold"  is  used  to  hold  the  value  of  the  previous  root  for  testing  against  the  present  ~  if 
N[Bold,3]  !=  N(B/.ruIe2,3]  the  root  is  considered  different  and  added  to  the  list  "Rulel" 
and  "rule2"  hold  the  results  of  "HndRoot"  for  the  three-and  four-layer  eigenvalue  equation, 
respecdvely.  Test  is  the  test  value  for  the  range  of  possible  roots.  List  is  the  intermediate 
list  of  output  values. 

3.1.4.  PROCEDURE  "PropagationConstant" 

"PropagadonConstant"  calculates  and  plots  all  confined  modes  for  the  input  parameters. 
For  each  propagadon  constant  found,  the  waveform  is  plotted  showing  the  mode  numbers 
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explicitly.  The  field  amplitude  coefficients  “ampTwoA,"  "ampTwoB,"  "ampThreeA," 
"ampThreeB,”  and  "ampFour."  are  calculated  by  Mathematica  using  the  eigenvalue 
equation  and  the  "Solve"  function;  "ampOne"  is  set  to  1.  The  function  "Equation"  is  the 
four-layer  equation,  k  is  the  propagation  constant,  increment  is  the  step  through  the  range 
of  possible  roots  and  list  is  the  list  of  propagation  constant  values.  The  parameters  q,  h,  p, 
and  s  are  the  same  as  those  defined  in  equation  (6). 
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PropagationConstant[l.,1.691, 1.46,3.85  -  0.02  1,1,1.43,1] 


10.0032  -  0.0000236124  I 
10.47  -  1.13667  10-6  I 

Begin  plotting  routine . 


Vave  Function 


Figure  2  Example  of  "PropagationConstant"  Procedure 
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3.1.5.  Attenuation  Calculation  Procedures 

"N1  Attenuation,"  "N2Attenuation,"  "N3Attenuation,"  "N4Attenuation," 
"GuideAttenuation,"  "CladdingAttenuation,"  and  "WavelcngthAttenuation"  all  have  the 
same  basic  outline.  Each  one  calculates  the  loss  due  to  attenuation  versus  a  varied 
waveguide  parameter.  For  a  description  of  the  procedure  see  "NlAttenuation." 


3.1.5.1  PROCEDURE  "  NlAttenuation" 

"NlAttenuation"  accepts  user  input  in  the  format: 

"NlAttenuation"!  { start,end,increment }  ,n2jn3,n4,cladding,guide,wavelength] 
where  {start,end,increment}  replaces  ni  and,  therefore,  defines  the  range  of  values  for  ni. 
All  equations  that  involve  ni  are  rewritten  replacing  ni  with  'variable'  where  'variable' 
ranges  fiom  start  to  end  in  steps  of  increment.  The  input  parameters  are  checked,  the  roots 
searched  for  and  the  plot  is  made.  All  parameters  and  local  variables  have  been  defined 
previously.  NumberQ[N[variab/€/uzme]]  is  used  to  check  the  input  before  any  processing 
has  taken  place  to  make  sure  the  variables  contain  numeric  values.  Since  a  user  enters  these 
values,  it  is  possible  that  the  values  are  not  numerical  as  they  should  be  and  any 
computation  will  crash. 
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NlAttenuanon[{l.,1.4.0.05},L67.1.46.3.85  -  0.02  1,0.5,0.5,0.6328] 


{1., -6.84822  10-*  } 
{1.05,-6.73274  10^  } 
(1.1,-6.60235  10 -«  } 
(1.15,  -6.45425  10^  } 
{1.2,  -6.28474  10  «  } 
{1.25,-6.08884  10^  } 
{1.3, -5.85966  10-«  } 
{1.35,-2.94539  10^  } 
{1.35,  -5.58731  10^  } 
{1.4,-2.39124  10^  ) 
{1.4,  -5.25691  10-*  } 


-8 

-5.25  10 

-8 

-5.75  10 

attenuation  (dB/cm)  vs  nl 

• 

• 

• 

-8 

1.1  1.2  •  1.3  1.4 

-6.25  10 

-8 

• 

-6.5  10  „ 

• 

-8 

• 

-6.75  10  . 

A.tteiiiiation  (dB/cm)  vs  xil 


Figure  3  Example  of  ”N1  Attenuation"  Procedure 
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3.1.5.2  PROCEDURE  "N2Attenuation” 


"N2Attenuation"  is  like  the  previous  attenuation  functions,  except  tha%  are  special 
cases  for  n2  since  it  may  be  a  complex  number.  The  variables  "imaginary Var"  and 
"realVar"  are  used  for  varying  such  parts  and  variable  is  the  combination  of  the  "Re"  and 
"Im"  part  A  variation  of  na  may  result  in  values  such  as  1.52  -  0.01 1  and  1.55  -  0.02  L 
The  attenuation  values  are  considered  different  for  the  number  of  modes,  but  the  values  will 


be  plotted  above/below  each  other  as  only  the  Reina]  is  used  in  plotting. 


N2Attenuation[l.,{1.6,1.68,0.02},1.46,3.85  -  0.02  1,0.5,0.5,0.6328] 

{1.6,-2.91513  10-''  } 

(1.62,  -1.90045  10-7  } 

(1.64,-1.25314  10-7  } 

{1.66,  -8.35437  10*8  } 

{1.68,  -5.6277  10*8  } 


Attenuation  (dB/cm^  vs  n2 


Figure  4  Example  of  "N2Attenuation"  Procedure 


3.I.5.3.  PROCEDURE  ”N3Attenuatioii" 


N3Attenuation[l.,1.67.{1.3,L5.0.05},3.85  -  0.02  1.0.5,0.5,0.6328] 

{1.3, -6.59959  10-<>  } 

(1.35,  -9.75759  lO*’  ) 

(1.35,-1.26601  10-«} 

(1.4,  -2.06399  10  } 

(1.4, -2.59717  10-«  } 

(1.45,-5.77298  10-«  } 

(1.5,  -1.41408  10-7  } 


Attexraation  (dB/cm)  vs  x3 


attenuation  (dB/cm)  vs  n3 


Figure  5  Example  of  "N3Attenuation"  Procedure 
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3.I.5.4.  PROCEDURE  “N4Attenuatioii" 


'*N4Attenuation”  is  like  the  previous  attenuation  functions,  except  there  are  special 
cases  for  n4  since  it  is  a  ctxnplex  number.  The  variables  "imaginaiy  Var"  and  "realVar"  are 
used  for  varying  such  parts  and  variable  is  the  combination  of  the  "Re"  and  "Im"  part  A 
variation  of  n4  may  result  in  values  such  as  3.80  -  0.01  I  and  3.80  -  0.02  I.  The 
attenuation  values  are  considered  different  for  the  number  of  modes,  but  the  values  will  be 
plotted  above/below  each  other  as  only  the  Re[n4]  is  used  in  plotting. 


N4Attenuation[l.,1.67,1.46,{3.80  -  0.01 1,3.9  -  0.02  1,0.05  -  0.005 
I), 0.5,0.5,0.6328] 


{3.8  -  0.01  1.  -6.95785  10-8  } 
(3.85  -  0.01  I,  -6.85647  10-8  j 
(3.9  -  0.01  1,  -6.75815  ICH  | 
{3.8  -  0.015  I,  -6.95353  10-8  } 
{3.85  -  0.015  I,  -6.85236  10-8  j 
{3.9  -  0.015  I,  -6.75423  10-8  } 
{3.8  -  0.02  I,  -6.94919  10-«  } 
{3.85  -  0.02  I,  -6.84822  10-«  } 
{3.9  -  0.02  I,  -6.75028  10-8  } 


-6.75  10 


-6.8  10 


-6.85  10 


-6.95  10 


Figure  6  Exan^le  of  "N4Attenuation"  Procedure 
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3.I.5.5.  PROCEDURE  "GuideAttenuation" 


GuideAttenuation[l.,1.67,1.46,3.85  -  0.02  1,1, {0.5,0.52,0.005}, 0.6328] 

(0.5,  -8.26635  10-“  } 

{0.505,-7.82954  10-“  } 

{0.51,  -7.42218  10-“  } 

{0.515,  -7.0419  10  “  1 
{0.52,  -6.43929  10-“  ) 


Attenuation  (dB/cm) 
vs 

Guiding  Layer  Thickness  (microns) 


Figure  7  Example  of  "GuideAttenuaticm''  Procedure 


21 


3.I.5.6.  PROCEDURE  "Cladding Attenuation" 


aaddmgAttenuation[1..1.67.1.46,3.85  -  0.02  I.{0.5,1.0.1 1,0.5,0.6328] 

{0.5,  -6.84822  10  «  ) 

(0.6,  -1.78486  10^  } 

(0.7,  -4.65518  10-«  } 

{0.8,  -1.21442  10-9  } 

{0.9,  -3.16838  10-10  } 


attenuation  (dB/cm) 

YS 

Cladding  Thickness  (microns) 


Figure  8  Exainpleof"CladdingAttenuation"  Procedure 
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3.I.5.7.  PROCEDURE  "WavelengthAttenuation" 


WaveleagthAttenuation[l.,1.67,1.46,3-85  -  0.02  I,0.5,0.5,{0.5,0.7,0.05}] 

{0.5,  -5.09033  10-''  } 

{0.55,  -1.70656  10^  } 

{0.6,  -4.13852  10-«  ) 

(0.65,  -8.73335  10-*  } 

(0.7,  -1.64989  10-7  } 


kttmuxmtxoxi.  (dB/cm) 

TS 

(nicro&s) 


Figure  9  Example  of  "WavelengtfaAttenuadon"  Procedure 
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3.1.6.  PROCEDURE  "CutoffThickness" 


The  cutoff  thickness  ("cutoffGuide")  is  calculated  from  the  three-layer  eigenvalue 
equation  as  the  value  of  the  effective  index  of  refraction  (proportional  to  the  propagation 
constant)  approaches  the  index  of  the  third  (cladding)  layer.  Error  messages  and 
appropriate  "Break"  statements  are  implemented  in  cases  of  error.  A  plot  of  cutoff 
thickness  versus  n2  is  output  along  with  the  values. 


CuloffThicknessK  l.S.1.6.0.01 )  .1..1.46j0.6328] 

1.5  0J68192 
IJl  0.319696 

1.52  0.28401S 

1.53  0.256386 

1.54  0.234201 
1J5  0.2159 

1.56  0.200485 

1.57  0.187282 

1.58  0.175819 

1.59  0.165753 

1.6  0.156829 


Cutoff  Thickness  (laicrons)  vs  n2 
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3.1.7.  PROCEDURE  "CutoffWavelength" 


The  cutoff  wavelength  (’’cutoffWave")  is  calculated  from  the  three-layer  eigenvalue 
equation  as  the  value  of  the  effective  index  of  refraction  approaches  the  index  of  the  third 
(cladding)  layer.  Error  messages  and  impropriate  "Break"  statements  are  inmlemented  in 
cases  of  oror.  A  plot  of  cutoff  wavelength  versus  n2  is  ou^ut  along  with  the  values. 


CutofrWavdength[(  l.S.l.6.0.01  }.1..1.46.0.S] 

1.5  0.859333 

1.51  0.98%92 

1.52  1.11402 

1.53  1.23408 

1.54  1.35098 

1.55  1.46549 

1.56  1J7817 

1.57  1.68943 

1.58  1.79958 

1.59  1.90886 

1.6  2.01748 


Cutoff  Vavelength  (laicrons)  vs  n2 


Figure  1 1  Example  of  "CutoffWavelength"  Procedure 
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3.2.  MODELING  EXAMPLE 

As  an  example,  consider  a  structure  consisting  of  0.4  pm  guiding  layer  of  PMMA 
(nsl.489)  on  0.65  pm  SiOa  (n=1.46)  on  a  silicon  substrate  (n=3.85-0.02 1)  operating  at  a 
wavelength  of  0.6328  pm.  It  is  useful  to  know,  before  actually  attempting  to  waveguide, 
how  many  guided  modes  are  supported  by  the  structure.  The  parameters  are  entered  into 
the  code  (see  Section  9.1.  for  details)  with  the  following  result: 

PropagationConstant[1.0,1.489,1.46,3.85  -  0.02  1,0.4,0.65,0.6328] 

There  are  no  guided  modes  for  this  four-layer  guide. 

The  structure  may,  however,  suppon  modes  at  other  wavelengths.  These  wavelengths  can 
be  located  using  this  code  as  outlined  in  Section  3.1.7.  The  result  is  shown  in  Hgure  12. 
It  can  be  seen  that  0.564pm  is  the  cucoff  wavelength  for  this  structure.  This  means  that  if  a 
shorter  wavelength  than  0.564pm  is  chosen,  there  is  a  guided  mode  for  this  structure.  We 
can  now  calculate  the  modal  behavior  of  the  guide  at  0.514pm.  The  propagation  constant  P 
and  the  fields  in  each  layer  are  calculated  and  ploned  as  shown  in  Figure  13.  Valuable 
experimentation  time  is  saved  by  noting  how  ill-couHned  this  particular  mode  is  and 
fabricate  another  sample  with  thicker  guiding  and  buffer  layers  for  better  ccmnnement  As 
a  final  step,  a  sample  with  0.5pm  PMMA  and  2pm  Si02  at  0.514pm  is  examined.  This 
calculation  is  shown  in  Hgure  14.  This  structure  has  two  guided  modes,  the  lower  order 
mode  being  fairly  well  confined. 
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CutoffWavelength[  { 1 .479, 1 .501,0.002 ) ,  1 .0, 1 .46,0.04] 


1.479  0.439212 
1.481  0.465665 
1.483  0.491297 
1.485  0.516219 
1.487  0.540522 
1.489  0.564279 
1.491  0.58755 
1.493  0.610388 
1.495  0.632833 
1.497  0.654925 
1.499  0.676694 


Cutoii  Vavelexigtlx  (microns)  vs  xk2 

0.6 
0.5 
0.4 
0.3 


1.46  1.47  1.48  1.49  Wavelength  (jim ) 

Guiding  Layer  Index  of  Refraction 

Figure  12  Theoretical  Calculation  of  Cutoff  Wavelength  for  0.4iim  PMMA  on 
0.65^  Si02  on  Si 
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PropagationConstant[1.0,1.489,L46,3.85  -  0.02  1,0.4,0.65,0.514] 
17.9227  -  0.00375371  1 
Begin  plotting  routine . 


V«.Te  Functioxi 


Distance  (pm) 


Figure  1 3  Theoretical  Calculation  of  the  Propagation  Constant  for  0.4pm  PMMA  on 
0.65pm  Si02  on  Si  at  2^.5 14mm 
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PropagationConstant[1.0,1.489,l-46,3.85  -  0.02  I,0.5,2,0.514j 

18.0052  -  0.000628331  1 
18.1522  -  0.0001144491 

Begin  plotting  routine . 


Amplitude  (Arbitrary  Units) 


Distance  (pm) 


Amplitude  (Arbitrary  Units) 

Figure  14  Theoretical  Calculation  of  the  Propagation  Constant  for  0.5pm  PMMA  on 
2.0pm  Si02  on  Si  at  X=0.514pm 
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4.0.  WAVEGUIDE  MATERIAL 


Recent  work  in  optical  waveguides  has  been  in  developing  a  material  that  has  both 
good  opdcal  qualities  as  well  as  long-tom  environmental  stability  and  durabilior.  One  area 
of  materials  interest  is  in  polymers.  Unfcntunately,  many  of  the  well-known  polymers  that 
exhibit  desirable  optical  characteristics  have  a  low  (less  than  100°C)  glass  transition 
temperature  (Tg)  which  can  cause  drastic  deterioration  of  several  key  mechanical  and 
physical  properties  as  the  polymer  chains  gain  high  mobility  at  moderate  temperatures. 
Low  Tg  polymers  are  also  more  susceptible  to  thermal  degradation  with  short,  high 
temperature  excursions  experienced  in  commercial  waveguide  processing  techniques. 
Having  a  high  Tgand  visual  transparency,  a  hexafluorc  isopropylidene-polybenzoxazole 
(6F-PBO),  became  a  primary  candidate. 

Having  been  synthesized  by  various  methods  approximately  25  years  ago^*^, 
polybenzoxazole  polymers  have  been  known  for  their  high  strength,  high  modulus,  and 
high  theimo-oxidadve  stability,  but  have  also  suffered  from  limited  solubility  and  low 
processability  which,  overall,  has  severely  hampered  their  material  development.  One 
manner  in  which  these  limitadons  have  been  overcome  is  through  the  incorporadon  of  a 
2,2-hexafluoroisopropylidene  (6F)  group  into  the  polybenzoxazole  backbone  to  produce  a 
more  transparent,  soluble  polymer  which  did  not  sacrifice  thermo-oxidadve  stability 
Having  high  use  temperature  comparable  to  that  of  epoxy  and  polyimide  systems,  the  6F 
polybenzoxazoles  (6F-PBO)  have  an  advantage  in  that  unlike  the  polyimides  which  are 
only  soluble  in  the  prepolymer  form,  it  is  soluble  in  common  organic  solvents  in  the  closed 
polymer  form.  They  also  have  the  advantage  that  their  Tg  can  be  varied  from  94®C  to 
382*’C  through  structural  modidcadon  of  the  aromadc  spacer  group  in  the  backbone^^.  The 
version  of  6F-PBO  which  was  examined  in  this  paper  contained  a  paralinked  diphenyl  ether 
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aromatic  spacer  group  (Figure  15)  which  resulted  in  a  Tg  of  300®C.  This  polymer  is 
denoted  6F-PBO-1E. 

Having  established  6F-PBO-1E  as  meeting  the  requirements  on  processability  and 
Tg,  the  optical  properties  were  then  examined.  This  technical  report  describes  our  findings 
in  the  investigation  of  the  optical  waveguide  loss  measurements  on  the  polymer  material 
6F-PBO-1E. 


PBO  (Poly-p-phenyleoe  benzobisthiazole) 
(rigid-r^  polymer) 


6F-PBO-1E  (End-capped) 


Figure  15  6F-PBO  and  Variations 

4.1.  WAVEGUIDE  FABRICATION 

4.1.1.  THE  SUBSTRATE 

The  substrate  material  used  for  the  optical  waveguide  fatmcadon  is  silicon.  Wafers 

with  <1(X)>  orientation  were  used.  Due  to  the  silicon  high  index  of  refraction,  an  oxide 

layer  is  necessary  for  guided  wave  propagation.  The  oxide  layer  was  grown  with  steam  at 
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1  lOO^C  for  approximately  4  hours.  Thickness  calculations  were  accomplished  with  the 
growth  modeling  package  TSupieme4  by  Technical  Modeling  Associates.  The  2-  and  3- 
inch  wafers  were  oxidized  separately  and  dieir  calculated  oxide  thicknesses  are  1.4  pm  and 
1.5  pm,  respectively. 

The  substrate  material  was  measured  for  thickness  and  refractive  index  values  using 
an  ellipsometry  analysis  software  package,  DAFEBM,  in  conjunction  with  a  fixed 
wavelength  (X  =  632.8nm)  Rudolph  Research  AutoEL  n  automated  nulling  ellipsometer. 
The  results  of  this  analysis  are  shown  in  Figures  16-19.  The  3-inch  wafers  have 
considerable  vaiiadon  in  both  index  and  thickness,  most  likely  due  to  a  slight  bowing  of  die 
wafer  after  oxidation.  Based  on  the  analysis  of  the  two  wafer  types,  the  majority  of  the 
guided  wave  loss  measurements  were  done  on  the  2-inch  wafers. 

Values  for  index  of  refraction  and  thickness  for  both  the  oxide  layer  and  the 
polymer  layer  are  averaged  over  the  2-inch  sample  values  only.  The  oxide  layer  thickness 
was  found  to  be  1.30pm  and  1.43pm  for  the  2-  and  3-inch  wafers,  respectively.  The 
index  of  refraction  of  the  oxide  layer  is  taken  to  be  1.455.  A  common  base  clean  procedure 
was  used  to  prepare  the  wafers  for  spinning  (see  Section  8.0). 


32 


Refractive  Index  g  Thickness  (A) 


Topical  3"  Wafa:  Oxide  Profile 


1 6  Oxide  Thickness  Profile  taken  on  Profilometer  for  3-Inch 
Oxidized  Silicon  Wafer 


Typical  3”  Wafer  Oxide  Profile 


Figure  17  Refractive  Index  Profile  of  Oxide  Layer  taken  on  Ellipsometer 
for  3-Inch  Oxidized  Silicon  Wafer 


index  G  Thickness  (A) 


Kgure  19 


Thickness  Profile  of  Oxide  on  2"  Wafer 


Oxide  Thickness  Profile  Taken  on  Profilraietry  for  2'Inch 
Oxidized  Silicon  Wafer 


Typical  Oxide  Index  Profile  on  2"  Wafers 


^  l-O  .JO 


0.0  - 

-0^ 


Reftactive  Index  Profile  of  Oxide  Layer  taken  on  Ellipsometa- 
for  2-Inch  Oxidized  Silicon  Wafer 
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4.1.2.  THE  GUIDING  MATERIAL 

The  polymer  was  prepared  by  Daychem,  Inc.  in  a  manner  similar  to  that  described 
by  Reinhardt^^  with  an  additional  post-polymerization  end-capping  of  the  polymer  chains 
with  benzoic  acid^^.  The  polymer  was  end-capped  to  inhibit  reaction  of  the  chain  ends  at 
elevated  temperatures.  Due  to  the  polymer  high  temperature  stability,  high  solubility,  and 
transparency,  it  appeared  to  offer  great  pronoise  as  an  optical  waveguide  material  and 
possibly  even  as  a  host  fcx*  a  guest-host,  second-order  nonlinear  optical  film. 

In  an  attempt  to  find  the  optimum  organic  solvent  that  would  both  dissolve  the 
polymer  and  produce  a  high  quality  spin-coated  film,  the  polymer  was  dissolved  in  several 
organic  solvents.  The  solvents  which  caused  the  most  solubility  included  chloroform, 
l,r,2,2'-tetrachloroethane.  2-picoline,  and  4-picoline  (each  at  5  wl%).  Each  solution  was 
pipetted  onto  1”  x  1"  microscope  glass  slide  and  spun  at  3(X)0rpm  for  60  seconds  on  the 
Solitec  5 1 10-CT  spin  coater  in  a  clean  room  environment.  Optical  inspection  and  Dektak  n 
2A  profilometry  results  indicate  that  the  film  spun  from  the  tetrachloroethane  solution  was 
the  most  transparent  film  with  the  iix)st  uniform  thickness  and  least  surface  roughness. 

Efforts  were  made  to  increase  the  concentration  of  the  tetrachloroethane  solution  in 
order  to  produce  a  thicker  film  under  the  same  spinning  conditions.  A  thicker  film  would 
support  more  modes  of  propagation  and  more  beam  confinement,  but  a  more  concentrated 
and  viscous  solution  would  also  be  more  difficult  to  syringe  filter  through  a  0.5|im  Teflon 
filter  as  is  required  for  reduction  of  dust  and  other  scattering  centers.  The  same  spin 
parameters  were  employed  due  to  the  initial  success  in  film  quality  and  thickness 
uniformity.  A  7  wt%  solution  of  the  6F-PBO-1E  in  tetrachloroethane  was  successfully 
filtered  after  many  unsuccessful  attempts  at  higher  concentrations. 


35 


The  6F-PB0-1E  from  the  7  wt%  solution  in  tetzachloroethane  was  spun  onto  the 
five  oxidized  silicon  substrates  at  3000  rpm  for  60  seconds.  Each  wafer  was  baked  out  at 
80°C  for  12  hours  following  the  spin  to  drive  off  the  tetrachloroethane  solvent 

4.2.  SAMPLE  CHARACTERIZATION 

Visual  inspection  of  the  resulting  baked  samples  revealed  limited  inclusions  and, 
overall,  relatively  good  film  quality.  These  samples  were  then  evaluated  to  determine  index 
of  refiacti<xi,  thickness,  and  absorption  before  attempting  to  waveguide. 

4.2.1.  GUIDING  LAYER  THICKNESS 

The  thickness  was  measured  on  a  Dektak  IIA  Profilometer.  A  radial  scratch  was 
induced  in  an  unaltered  region  of  the  polymer  film  and  the  step-height  of  the  scratch  was 
measured.  Figure  20  shows  the  results  of  this  study.  The  average  over  all  five  san^les 
was  1.00^m  and  this  value  was  used  in  the  theoretical  evaluation.  As  is  typical  in  spin 
coating  procedures,  there  is  a  relatively  ui^orm  thin  area  in  the  center  of  the  sample,  and 
increasing  thickness  at  the  edges.  Due  to  this  fact,  all  measurements  were  made  on  the 
uniform  portion  of  the  samples.  In  addition,  this  flat  uniform  region  is  more  prominent  on 
the  2-inch  wafers,  and  the  infotination  previously  presented  for  the  oxide  layer  and  index 
of  refiraction  data  was  only  recorded  on  those  samples. 
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11000 


Typical  Thiclaiess  Profile 
(Scan  on  2”  wafer  by  Profilometer) 


Hgure20  Polymer  Thickness  Rofiletaten  on  Plrofilomctcr  for  2-Inch  Wafer 


4.2.2.  INDEX  OF  REFRACTION 

The  index  of  refraction  of  the  polymer  film  was  determined  using  a  four-layer 
(substrate  -  oxide  -  film  -  air)  ellipsometry  analysis  software  package,  DAFIBM,  in 
conjunction  with  a  fixed  wavelength  (X  =  632.8nm)  Rudolph  Research  AutoEL  II 
automated  nulling  ellipsometer.  If  the  index  and  thickness  of  the  oxide  layer  and  the 
thickness  of  the  polymer  layer  is  known,  this  software  program  can  be  used  to  calculate  the 
index  of  the  polymer  layer. 

The  index  of  refraction  was  found  by  performing  two  perpendicular  radial  scans  on 
both  the  2-  andS-inch  wafers.  Figure  21  shows  the  data  taken  on  the  2-inch  sample  with 
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the  ellipsometer.  Analysis  of  the  data  showed  that  a  variation  in  index  would  generally 
consist  of  a  level  region  near  1.69,  increasing  to  indices  close  to  2.0  on  the  edges.  The 
higher  index  outo:  regions  are  currently  a  source  of  further  study. 

To  ctmfinn  the  index  measurements  made  using  the  ellipsometer,  a  more  easily 
analyzed  single  layer  of  polymer  on  bare  silicon  was  fabricated  under  the  same  processing 
conditions  from  the  samg  solutimi  as  the  waveguiding  samples.  In  addition  to  analysis  by 
ellipsometry,  this  new  film  was  analyzed  using  a  Metricon  2010  Prism  Cmipler  where  the 
index  and  thickness  of  the  film  were  measured  by  recording  the  angles  at  which  modes  are 
coupled  into  and  out  of  the  film  by  a  synunetric  prism  (technique  described  by  Tien  and 
Ulrich^^),  Index  and  thickness  measurements  were  obtained  in  the  iimer  1-inch  diameter  of 
a  2-inch  sample.  Within  this  region,  the  average  index  as  measured  with  the  ellipsometer 
was  found  to  be  1.668  (2.5%  variation),  while  over  five  equally  spaced  regions  on  the 
wafer,  the  Metricon  measured  the  index  to  be  1.6647  (0.01%  variation)  in  the  TE  irxxle  at 
632.8nm. 

The  average  polymer  layer  index  of  refraction,  over  the  three  procedures  perfcnmed 
(see  Table  1),  is  1.691. 


38 


Refractive  Index  Profile  in  Uniform  Region  of  2"  Wafer 
(Scanned  by  Ellipsometer) 


Hguie  21  Refractive  Index  Profile  of  Polymer  Layer  taken  on  Ellipsometer 
for  2-Inch  Wafer 


Procedure 

«F-PBO-ECoa 
Oxidized  Wafers 

6F-PBO-ECon 
UiMMddized  Wafers 

Qlipsoaeicr 

Index 

TWckaai 

Index 

TMckacv 

1^ 

aSSilia 

I.6a 

l.06)im 

Menieoa 

1.6647 

I.06|tm 

— 

LOOani 

— 

TABLE  1  Comparison  of  Three  Methods  for  Finding  Index  of  Refraction 
and  Thickness  of  6F-PBO-1E 
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4.23.  ABSORPTION 


Material-inherent  properties  are  not  considered  in  the  theoretical  calculations  so  it  is 
appropriate  to  investigate  the  effects  the  material  itself  will  have  on  the  waveguiding 
properties. 

In  order  to  evaluate  the  absorption  losses  for  6F-PBO-1E,  a  free-standing  film  was 
cast  through  slow  evaporation  of  the  solution  in  a  Petri  dish.  Transmission  studies  of  this 
film  at  normal  incidence  over  the  wavelength  range  250  -  320Qnm  were  then  performed 
using  a  Hitachi  U-4001  UV/VIS/NIR  spectrophotometer.  The  results  are  shown  in  Hgure 
22.  For  an  absorbing  film  at  ntmnal  incidence,  the  following  equations  for  Reflectivity, 
R,  Transmissivity,  T,  and  absorption  coefficient,  a,  apply.^^’*^ 

fl  -  RV^  e-«<* 

T  =  ^ -  where  d  is  the  thickness  of  the  film  (19) 

1  -  R2e-2a<» 

a-ll^(l-RF-^V(l^-Rf +  4T^r^|  (20) 

Assuming  the  film  to  be  non-absorbing  (a  =  0)  in  the  region  where  the  transmission  is 
maximum,  the  index  of  the  film  can  be  calculated  using^^= 

T  =  (21) 

1  -  R2 


where  the  reflectance  for  a  free-standing,  non-absorbing  film  is^= 


Transmission  Spectra  for  Free-Standing  Film 


Figure  22  Transmission  Data  Taken  on  Free-Standing  Film  of  6F-PBO  IE 
with  Spectrophotometer 


Solving  for  the  index  of  refraction  gives; 

2  -4.  Y4  -  4  T2  (23) 

2T 

The  assumption  is  then  made  that  dispersion  is  negligible  so  the  index  found  at  the 
maximum  transmission  wavelength  is  a  constant  in  the  calculation  of  the  absorption 
coefficient  over  all  wavelengths.  The  reflectivity  is  also  assumed  to  be  a  constant.  These 
assumptions  allow  only  for  a  semi-qualitadve  analysis  of  the  film,  but  this  is  sufficient  to 
investigate  the  absorption  trends. 
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The  absorption  coefficient  a  can  then  be  calculated  using: 

a  .  1.  (24) 

where  R  utilizes  the  formula  for  non-absorbing  films,  T  is  that  measured  experimentally 
with  the  spectrophotometer,  and  d  is  the  thickness  of  the  film  as  measured  with  the 
profilometer. 

Figure  23  depicts  the  findings  of  this  investigation.  The  loss  coefficient,  a,  as 
calculated  from  equation  (23)  depicts  strong  absorption  in  the  shorter  wavelength  region. 
This  information  was  useful  in  evaluating  the  experimental  waveguide  loss  data. 

Absorption  Coefficient  of  Free-Standing  6F-PBO-1E  Film 


Figure  23  Absorption  (Coefficient  for  6F-PBO-1E  as  Calculated  from 
the  Transmission  Data 
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4.3.  THEORETICAL  BEHAVIOR 


Once  the  index  of  refraction  and  thickness  of  both  the  oxide  and  guiding  layers  of 
the  samples  have  been  measured,  the  package  presented  in  Section  3.0  was  used  to  predict 
the  behavior  of  the  material  as  a  planar  waveguide. 

The  loss  of  each  confined  mode  has  been  calculated  at  several  wavelengths,  and  the 
combined  results  ere  shown  in  Figure  24.  There  is  at  least  one  confined  mode  for  all 
wavelengths  of  interest,  therefore,  the  samples  fabricated  theoretically  fulfill  the 
requirements  for  measuring  the  opdeal  attenuation. 


Theoretical  Loss  versus  Wavelength 
I.OOiim  6F-PBO-1E  on  1.43pm  SIO^  on  SI 


Hguie  24  Theoretical  Plot  of  Loss  versus  Wavelength  Computed  with  Maihematica 
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5.0.  EXPERIMENTAL  MEASUREMENTS  OF  WAVEGUIDE 

ATTENUATION 


5.1.  Introduction 


When  studying  a  material  for  possible  use  in  optical  integrated  circuits,  one 
important  consideration  is  attenuation  in  the  waveguide.  Waveguide  attenuation  may  have 
many  causes  including  absorption,  scattering,  and  substrate  coupling. 

The  energy  loss  in  the  waveguide  may  be  modeled  using  Beer's  Law; 


I  =  lo  e*®  ^  where  I©  is  the  input  intensity  (25) 

a  is  the  loss  coefficient  and 
z  is  the  distance  of  propagation 


By  definition: 

IdBslOlogj^J 


(26) 


Taking  z  =  1  cm  and  transforming  logarithm  bases  yields: 
log  =  -0.434  a  cm 


(27) 


Transforming  this  expression  into  dB  notation  gives: 
#  =  -4.34  a 


where  a  is  die  attenuation  coefficient. 


(28) 
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5.2.  LOSS  MEASUREMENT  TECHNIQUES 

The  absoiption  coefficient,  a,  can  be  measured  by  any  of  three  common  loss 
measurement  methods:  three-prism  method^^,  two-prism  method,  or  out-of-plane 
scattering  method  This  study  uses  all  three  and  compares  the  results.  A  custom-designed 
waveguide-loss-measurement  apparatus^^  is  used  to  perform  all  three  types  of 
measurements. 

5.2.1.  THE  THREE  PRISM  METHOD 

In  the  three  prism  method^^,  the  input,  coupling,  and  output  prisms  and  their 
associated  detectors  are  arranged  as  in  Figure  25.  Prisms  1  and  3  are  clamped  while  2  is 
moveable.  Light  is  coupled  into  the  guide  with  prism  1,  called  the  input  prism,  and  out  of 


Prism  #2 

Measurement  Prism 


Prism  #3 
Exit  Prism 


Hgure  25  Experimental  Setup  for  the  Three-Prism  Method  to  Measure  Attenuation 
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the  guide  through  prisms  2  and  3,  referred  to  as  the  measurement  and  output  prisms, 
respectively.  Measurement  of  the  light  intensity  from  both  prisms  2  and  3  eliminates  the 
need  to  know  the  coupling  coefficient  for  the  moveable  piism^^. 

The  setup  works  by  first  measuring  the  output  at  prism  3  without  prism  2  clamped,  then 
measuring  the  output  of  the  2nd  and  3rd  prisms  as  a  function  of  distance  as  prism  2  moves 
along  the  guide.  The  intensity  can  then  be  calculated  as  a  function  of  distance  with  the 
relation: 


p9-P3 


(29) 


where  P3O  is  the  output  measured  at  prism  3  without  prism  2  clamped,  P2(2),  and  P3  are 
the  outputs  of  prisms  2  and  3,  respectively. 

A  Schott  SF-6  45“-45®-90®  coupling  prism  and  a  Hamamatsu  S1227-66BQ 
photodiode  with  a  33  mm^  active  area  are  used  in  each  of  the  measurement  positions.  The 
active  area  of  the  detector  is  parallel  to,  and  covers  the  out-coupling  face  of,  the  prism.  The 
prism-detector  assembly  moves  in  unison  to  insist  on  uniformity.  The  procedure  is 
computer  controlled  using  a  Hewlett-Packard  9000-330  computer  and  multiprogrammcr. 
The  waveguide  is  placed  on  the  platform,  and  the  input  coupling  prism  is  clamped  in  place. 
The  operator  adjusts  the  input  beam  for  propagation  through,  and  total  internal  reflection 
(TTR)  inside,  the  coupling  prism.  TTR,  by  definition,  creates  evanescent  fields  outside  the 
prism  boundaries,  in  the  air  gap  between  the  prism  and  sample.  This  evanescent  freld  is 
coupled  into  the  sample  and  becomes  a  guided  mode  if  TTR  is  also  satisfied  inside  the 
polymer  layer.  The  measurement  and  ou^ut  coupling  prism  assemblies  are  then  positioned 
over  the  guide.  The  computer  records  the  clamping  pressure  by  means  of  load  cells  located 
in  the  prism  arm  assembly.  The  operator  chooses  a  measurement  distance  and  increment 
size  over  which  to  take  measurements  and  then,  subsequently,  clamps  and  unclamps  the 
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measuiement  prism  as  the  computer  moves  the  measurement  arm  and  records  the  data. 
Finally,  measured  intensity  data  are  output  m  a  ccxnputer  file  for  plotting. 

Once  the  intensity  profile  has  been  established,  the  data  can  be  Htted  to  an 
exponential  function  (according  to  Beer's  Law),  and  the  attenuation  coefficient  can  be 
determined. 

5.2.2.  THE  TWO  PRISM  METHOD 

The  two-piism  loss  measurements  are  performed  in  much  the  same  manner  as  the 
three-prism  method,  except  of  course,  the  elimination  of  the  need  for  the  third  prism.  In 
the  two  prism  method,  the  load  cells  incorpcxated  in  the  prism  arm  assembly  of  the  custom 
apparatus^^  allow  assurance  of  constant  coupling  efficiency,  therefore,  eliminating  the  need 
for  the  third  prism.  One  disadvantage  in  both  the  prism  methods,  however,  com  ■  when 
measuring  soft  films.  Each  time  the  prism  is  clamped  down  on  the  sample,  it  has  the 
possibility  of  damaging  the  Elm  surface.  For  this  reason,  a  non-contact  method  of  loss 
measurement  is  primarily  used  in  this  study. 

5.2.3.  THE  OUT-OF-PLANE  SCATTERING  METHOD 

The  non-contact  method  for  determining  the  attenuation  coefficient  is  the  third 

cotrunon  waveguide  loss  measurement  technique,  referred  to  as  the  out-of-plane  scattering 

method.  In  this  method,  there  is  an  input  coupling  prism  and  a  moveable  detector  as 

shown  in  Figure  26.  The  detector  is  a  Hamamatsu  S1227-66BQ  photodiode  with  .a  33 

mnP  active  area  masked  to  a  1 -mm- wide  slit  positioned  perpendicular  to  the  guiding  streak. 

The  detector  is  moved  along  the  guiding  streak,  holding  the  position  above  the  sample 

constant  with  the  micrometer  and  arm  assembly.  The  intensity  of  the  light  scattered 

perpendicular  to  the  guiding  streak  is  measured  as  a  function  of  distance  from  the  coupling 

47 


spot  The  intensity  profile  is  fit  to  an  exponential  function,  and  the  attenuadtxi  coefficient  is 
determined  from  the  slope  of  the  curve.  This  method  is  much  simpler  than  the  three-piism 
method,  and  is  not  only  used  with  delicate  films,  but  with  all  types  of  samples. 

The  apparatus  used  for  the  out-of-plane  scattering  method  is  the  same  as  that  used 
for  the  prism  methods.  In  this  procedure,  the  first  two  arms  are  used.  The  measurement 
arm  no  longer  has  a  prism  attached  but  is  fit  with  a  detector  which  is  parallel  to  the 
platform.  As  in  the  three-prism  method,  the  operator  adjusts  the  input  beam  for  TIR. 
However,  the  computer  now  controls  the  movement  of  the  detector  along  the  streak, 
records  all  pertinent  information,  and  outputs  it  to  a  file  for  inspection. 


Hguie  26  Experimental  Setup  for  the  Out-of-Plane  Scattering  Method, 
Showing  Detector  Positioning  and  Masking. 
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5.3.  EXPERIMENTAL  PARAMETERS 


Loss  measurements  were  conducted  at  wavelengths  between  0.4S79  and  1.0|im. 
The  sources  were  a  Spectra-Physics  Argon  Ion  laser,  a  Coherent  599-01  Dye  laser  using 
Rhodamine  6G  and  a  Spectra-Physics  3950  TirSapphire  laser.  A  Burleigh  Model  WA-10 
Wavemeter  WA-10  was  used  to  measure  the  wavelength  of  the  input  light  and  a  Newport 
SP  550  polarization  rotator  was  used  to  ensure  TE  modes.  Schott  SF-6  glass  (n= 1.805) 
prisms  were  used  to  couple  the  light  into  the  guide.  The  beam  was  weakly  focused  at  the 
point  of  contact  with  the  guide  with  a  1000  mm  focal  length  lens.  Loss  measurements  were 
performed  on  each  of  five  samples  over  the  wavelength  range  0.46  -  0.975|im. 

Five  samples  of  end-capped  6F-PBO  were  fabricated  on  oxidized  silicon  wafers,  two  on 
3-inch  wafers  and  three  on  2-inch  wafers.  Preliminary  observations  of  guided  streaks 
suggest  this  material  has  a  relatively  high  loss.  Therefore,  most  of  the  loss  measurements 
were  made  using  the  out-of-plane  scattering  techiuque  on  the  2-inch  wafers.  However, 
some  three-prism  data  were  taken  on  the  3-inch  wafers.  The  2-inch  wafers  could  not  be 
measured  with  the  three-prism  method  due  to  a  size  ccmstraint  Some  data  were  taken  on 
these  samples  by  the  two-prism  method  where  the  coupling  coefficient  of  tiie  measurement 
prism  was  strictly  monitored  and  kept  constant  throughout  the  experiment  using  the  load 
cells  accommodated  in  the  arm  of  the  stage.  Due  to  the  destructive  nature  of  the  prism  out- 
coupling  methods  on  polymer  materials,  very  few  data  points  could  be  collected  on  these 
sair^les  and  only  at  wavelengths  where  the  material  has  low  absorption. 
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5.4.  RAW  WAVEGUIDE  DATA 

Hguies  27-32  are  representative  plots  of  the  raw  data  taken  on  one  of  the  2-inch 
samples  with  the  out-of-plane  scattering  method.  The  raw  data  consists  of  the  relative 
positions  of  the  measurement  prism  and  die  corresponding  intensities.  The  data  are  plotted 
with  the  software  package  KaleidaCraph  Version  2.1.3,  and  the  standard  exponoitial  fit 
included  in  that  package  is  calculated.  The  best  fit  exponential  curve,  along  with  the 
coireladon  value,  R,  is  displayed  in  the  bottom  left  corner  of  each  plot.  The  san^ile  used  to 
obtain  the  data  in  Figures  27-32  consists  of  a  LOpm  ftlm  of  6F-PBO-1E  on  1.4p.m 
oxidized  silicon  on  a  2-inch  silicon  wafer.  The  operating  wavelengths  are  noted  (m  the 
plots  along  with  the  measured  loss  in  dB/cm. 


6F-PBO-1E  on  Oxidizod  Silicon 


Rgure  27  Raw  Data  for  1 .0pm  6F-PBO-  IE  on  1 .4fim  Si02  on  Si  at  488nm 
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6F-PB0>1E  on  Oxidizod  Silicon 


nMBM  rVBIBm  (HIHIp 

Hgure  28  Raw  Data  for  l.Oum  6F-PB0-1E  on  1.4mn  SiC)2  on  Si  at  514.5nm 


6F>PBO>1E  on  Oxidizod  Silicon 


Figure  29  Raw  Data  for  l.O^m  6F-PBO-1E  on  1.4^m  Si02  on  Si  at  578.15nm 


6F>PBO-1E  on  Oxidizod  Silicon 


Figure  30  Raw  Data  for  l.O^m  6F-PB01E  on  1.4^  Si02  on  Si  at  SOOnm 


6F>PBO-1E  on  Oxidizod  Silicon 


Figure  31  Raw  Data  for  l.O^m  6F-PBO-1E  on  1.4ixm  Si02  on  Si  at  900nm 


6F-PB0<1E  on  Oxidizod  Silicon 


Hguxe  32  Raw  Data  for  l.CHtm  SF-PBO-IE  on  1.4tun  SiC)2  on  Si  at  975nm 


5.5.  THREE  METHODS  COMPARISON 

Figure  33  shows  the  comparison  on  the  three  methods  of  loss  measurement  on  the  3- 
inch  wafers.  The  sample  used  to  obtain  this  data  consists  of  a  1.0^m  film  6F-PBO-1E  on 
1.5^m  oxidized  silicon  on  a  3-inch  silicon  wafer.  The  abs(»ptic»i  coefScient  of 
the  3-inch  wafers  are  somewhat  lower  than  that  of  the  2-inch  wafers.  This  is  true  of  all 
samples  tested,  probably  due  to  the  substrate  irregularity  as  pointed  out  in  Section  4.1.1. 
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Comparison  of  Three  Methods  of  Loss  Measurements 
6F-PBO-1E  on  Silicon  Wafers 


Figure  33  Comparison  of  Three  Methods  of  Loss  Measurement  on  1 .0pm 
6F-PBO-1E  on  1.4pm  Si02  on  3"  Si  Wafer  at  900nm 


5.6.  WAVEGUIDE  RESULTS 

The  average  values  for  optical  waveguide  attenuation  in  dB/cm  along  with  their 
standard  deviations  are  shown  in  Rgure  34  and  TABLE  2. 

Upon  inspection  of  the  experimental  data,  it  is  observed  that  the  loss  increases  at 
both  the  long  and  shon  wavelength  ends  of  the  tested  spectrum.  The  increase  at  the  long 
wavelength  end  agrees  with  what  would  be  expected  since,  as  the  wavelength  increases, 
the  confinement  decreases  which  results  in  more  coupling  into  the  buffer  and  substrate 
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layers  and  mcne  energy  lost  to  the  system.  The  high  losses  at  shorter  wavelengths  may  be 
due  to  scattering  and  a  significant  increase  in  absoipdon.  Scattering  losses  are  expected  to 
increase  as  the  wavelength  shortens  due  to  the  fact  that  the  wavelength  is  approaching  the 
particle  size.  In  addition,  the  material  itself  is  more  absorbing  at  the  lower  wavelengths. 
This  agrees  with  what  was  determined  experimentally. 


Loss  versus  Wavelength* 
6F-PBO-1E  on  Silicon  Wafer 


0.4  0.5  0.6  0.7  0.8  0.9  1 

Wav«(6figlh  (microns) 

'Sunmiaiy  d  all  dau  points  mrilh  STD<0.5 


Figure  34  Summary  of  Loss  Data  for  6F-PBO-1E  at  0.45  -  0.975^m  Including 
Out-of-Plane,  Three-Prism  and  Two-Prism  Methods 
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2* 

AlloMthods 

3*  SsunplfiS. 
AUnelhods 

S'Sampia 

3  Priim  Method 

Wavelength  (nm) 

Lon 

inp 

Lou 

MMW 

*573 

0.61 

B 

H 

1.13 

0566 

B 

B 

AUJO 

5^40 

2505 

— 

B 

OM 

— 

— 

S01.7 

654 

aii3 

— 

— 

— 

SU3 

7J2 

ai36 

— 

— 

157 

— 

72S.0 

4.19 

1.137 

— 

— 

6.89 

7565 

— 

— 

7S0j0 

4.60 

1.476 

— 

— 

2.17 

— 

— 

— 

775.0 

SJO 

0.430 

— 

— 

— 

— 

— 

800i> 

457 

0531 

— 

B 

B 

— 

— 

82S.0 

5.71 

0.404 

— 

— 

— 

— 

— 

aSOjO 

5.71 

2592 

553 

RSI 

6.62 

2598 

— 

— 

S75.0 

756 

^545 

in 

mi 

654 

1588 

— 

— 

9oao 

6Jt4 

2556 

4.12 

0599 

5.47 

1.646 

4.38 

0539 

925.0 

757 

EE31 

— 

— 

5.18 

1.634 

4.84 

Wf!l 

9sao 

752 

0598 

— 

555 

1.460 

wm 

0.252 

975.0 

757 

0598 

B 

B 

650 

Z1S7 

3.79 

0.851 

TABLE  2  Summary  of  Average  Loss  Data  for  6F-PBO-1E 
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6.0.  CONCLUSION 


A  modeling  package  that  can  be  used  in  conjunction  with  material  fabricadons  to 
evaluate  possible  new  optical  waveguide  materials  has  been  discussed.  In  addition,  the 
polymer,  6F-PBO-1E,  has  been  investigated  to  ascertain  its  eligibility  as  a  waveguiding 
material.  The  subsequent  results  of  the  characterization,  modeling,  and  experimentation  of 
6F-PB01E  after  fabrication  into  waveguide  form  have  been  reported. 

The  polymer  6F-PBO-1E  has  been  characterized  to  determine  its  index  of 
refraction,  absorption,  and  optical  waveguide  attenuation.  The  measured  index  of  1.691, 
along  with  its  low  absorption  at  longer  wavelengths,  nuit.es  6F-PBO-1E  attractive  as  an 
optical  waveguide  material.  However,  with  losses  ranging  from  4  to  9  dB/cm  over  the 
visible  wavelengths,  while  it  may  still  be  useful  in  other  applications,  this  material  does  not 
meet  the  Air  Force  (e.g.  WL/ELOT,  RL/ERO)  loss  guideline  of  1  dB/cm  or  less  fw  optical 
applications.  We  conclude  that  this  material,  6F-PBO>lE,  in  its  present  form,  does  not 
warrant  further  investigation  as  an  Air  Force  applicable  optical  waveguide  material. 

Current  work  with  the  hexafluoro-isopropylidene-polybenzoxazole  polymers, 
however,  involves  the  investigation  of  variations  on  the  6F-PBO-1E,  via  changing  the  end¬ 
capping  material,  as  a  viable  candidate  for  high  temperature  second-order  nonlinear  optical 
host  material. 
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APPENDIX  A:  BASE  CLEAN  PROCEDURE 


20 


Materials:  Beaker  of  the  appn^ate  size 

Ammonium  Hydroxide  (NH4OH) 

Deionized  water 

Hydrogen  Peroxide  (30%) 

Thermometer 

Hotplate 


1 .  Carefully  mix  a  4: 1  H20:NH40H  solution  in  the  beaker. 

Add  NH4OH  to  the  H2O  slowly. 

2.  Heat  mixture  to  75®C. 

3 .  Add  the  same  quantity  of  hydrogen  peroxide  as  NH4OH 
to  the  solution. 

4.  Slowly  dip  the  wafer  in  the  beaker  and  allow  it  to  soak  in  the 
solution  for  IS  minutes,  maintaining  the  temperature  at 
65-75°C. 

5 .  At  the  end  of  IS  minutes,  put  the  entire  beaker  under  a  running 
stream  of  deionized  water  for  S  minutes. 

6.  Remove  the  wafer  and  blow  dry  with  dry  Nitrogen. 
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APPENDIX  B:  RAW  WAVEGUIDE  DATA  TABLES 


Out-of-plane  scattering  method,  first  3"  wafer 


Increment 

(mm) 

Total  Distance 
(mini 

Wavelength 

fnml 

-Alpha 

R-value 

0.10 

2 

457.9 

0.10 

2 

457.9 

0.0374 

0.980 

0.10 

2 

457.9 

0.0245 

0.950 

0.10 

2 

472.7 

0.0286 

0.905 

0.10 

2 

472.7 

0.0135 

0.892 

0.10 

2 

476.5 

0.0112 

0.887 

0.10 

2 

488.0 

0.0193 

0.912 

0.10 

1 

496.5 

0.0365 

0.883 

0.10 

1 

496.5 

0.0336 

0.936 

0.10 

1 

514.5 

0.0292 

0.942 

0.10 

1 

528.7 

0.0345 

0.941 

0.10 

2 

528.7 

0.1912 

0.917 

0.10 

2 

702.5 

0.1924 

0.999 

0.20 

3 

725.0 

0.4264 

0.912 

0.20 

3 

725.0 

0.2204 

0.894 

0.10 

4 

692.5 

0.0157 

0.918 

0.10 

2 

692.5 

0.0244 

0.869 

0.10 

3 

698.0 

0.0428 

0.983 

0.10 

4 

698.0 

0.0421 

0.967 

0.10 

4 

698.0 

0.0464 

0.962 

0.10 

4 

696.0 

0.0425 

0.965 

0.10 

4 

696.0 

0.0557 

0.841 

0.10 

4 

696.0 

0.0394 

0.923 

0.10 

3 

700.0 

0.0883 

0.922 

0.10 

3 

700.0 

0.0629 

0.940 

0.10 

3 

725.0 

0.0415 

0.948 

0.10 

3 

725.0 

0.0398 

0.960 

0.10 

3 

725.0 

0.0654 

0.980 

0.10 

3 

750.0 

0.0501 

0.981 

Out-of-plane  scattering  method,  first  2"  wafer 


Increment 

./mm) 

Total  Distance 
fmm) 

Wavelength 

(nm) 

-Alpha 

R-value 

0.50 

25 

705.0 

0.0217 

0.993 

0.50 

25 

705.0 

0.0225 

0.989 

0.50 

25 

725.0 

0.0260 

0.993 
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Inorement 

(mm^ 

Total  Distance 
(mml 

Wavelength 

fnm^ 

-Alpha 

R-value 

0.25 

15 

750.0 

0.0476 

0.997 

0.25 

15 

750.0 

0.0409 

0.993 

0.20 

10 

718.8 

0.0905 

0.988 

0.20 

10 

718.8 

0.1121 

0.992 

0.20 

10 

710.0 

0.1166 

0.991 

0.20 

10 

725.0 

0.1118 

0.992 

0.20 

10 

725.0 

0.1077 

0.994 

0.20 

10 

725.0 

0.1111 

0.993 

0.20 

10 

733.0 

0.1143 

0.988 

0.20 

10 

733.0 

0.1172 

0.987 

0.20 

10 

733.0 

0.1150 

0.987 

0.20 

10 

750.0 

0.1283 

0.991 

0.25 

10 

750.0 

0.1280 

0.992 

0.25 

10 

750.0 

0.1376 

0.994 

0.25 

10 

765.0 

0.1364 

0.992 

0.25 

10 

765.0 

0.1438 

0.987 

0.25 

10 

775.0 

0.1327 

0.990 

0.25 

10 

775.0 

0.1371 

0.992 

0.25 

10 

785.0 

0.1323 

0.990 

025 

10 

785.0 

0.1296 

0.989 

0.25 

10 

800.0 

0.1209 

0.975 

0.25 

10 

800.0 

0.1229 

0.969 

0.25 

10 

710.0 

0.0937 

0.989 

0.25 

10 

725.0 

0.0913 

0.989 

0.25 

10 

733.0 

0.1047 

0.981 

0.25 

10 

750.0 

0.1215 

0.970 

0.25 

10 

765.0 

0.1266 

0.976 

0.25 

10 

775.0 

0.1257 

0.979 

0.25 

10 

785.0 

0.0805 

0.974 

0.25 

10 

710.0 

0.1271 

0.971 

0.20 

10 

710.0 

0.1093 

0.993 

0.20 

10 

710.0 

0.1053 

0.983 

0.20 

10 

710.0 

0.1073 

0.985 

0.20 

10 

765.0 

0.1341 

0.992 

0.20 

10 

775.0 

0.1300 

0.990 

0.20 

10 

785.0 

0.1267 

0.986 

0.20 

10 

800.0 

0.1081 

0.986 

0.20 

10 

815.0 

0.1070 

0.987 

0.20 

10 

815.0 

0.1045 

0.987 

0.20 

10 

815.0 

0.1188 

0.962 

0.10 

8 

825.0 

0.1314 

0.973 

0.15 

8 

825.0 

0.1334 

0.974 

0.15 

7 

825.0 

0.1461 

0.980 

0.15 

8 

833.0 

0.1444 

0.978 

0.15 

8 

833.0 

0.1441 

0.978 

0.15 

8 

833.0 

0.1421 

0.975 

0.15 

8 

850.0 

0.1325 

0.981 

0.15 

8 

850.0 

0.1316 

0.979 

0.15 

8 

850.0 

0.1398 

0.982 
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Out-of'plane  scattering  method.  Second  2"  wafer 


Increment 

fmm) 

Total  Distance 
(mm) 

Wavelength 

(nml 

-Alpha 

0.25 

10 

710.0 

0.0878 

0.25 

10 

710.0 

0.0914 

0.15 

8 

725.0 

0.1123 

0.15 

8 

725.0 

0.1150 

0.15 

8 

733.0 

0.1148 

0.15 

8 

733.0 

0.1121 

0.25 

10 

750.0 

0.1048 

0.25 

10 

765.0 

0.1006 

0.25 

10 

765.0 

0.1102 

0.25 

10 

775.0 

0.1062 

0.25 

10 

775.0 

0.1168 

0.25 

10 

785.0 

0.1079 

0.25 

10 

785.0 

0.1109 

0.15 

8 

800.0 

0.1221 

0.15 

8 

815.0 

0.1177 

15.00 

8 

815.0 

0.1162 

0.15 

8 

800.0 

0.1209 

0.15 

8 

825.0 

0.1231 

0.15 

8 

825.0 

0.1239 

0.20 

8 

833.0 

0.1301 

0.25 

8 

833.0 

0.1270 

0.25 

8 

850.0 

0.1374 

0.25 

8 

850.0 

0.1330 

0.10 

5 

572.3 

0.0489 

0.10 

5 

572.3 

0.0491 

0.10 

5 

572.3 

0.0738 

0.10 

5 

578.2 

0.1050 

0.10 

5 

578.2 

0.0951 

0.15 

7 

578.2 

0.0497 

0.15 

7 

584.7 

0.0999 

0.15 

7 

584.7 

0.9061 

0.15 

7 

584.7 

0.0917 

0.15 

7 

590.0 

0.0949 

0.15 

7 

590.0 

0.1006 

0.10 

5 

590.0 

0.1089 

0.10 

5 

593.9 

0.1045 

0.10 

5 

593.9 

0.0927 

0.10 

5 

593.9 

0.1158 

0.10 

4 

528.7 

0.1545 

0.10 

4 

582.7 

0.1456 

0.10 

4 

514.5 

0.1702 

0.10 

4 

514.5 

0.1708 

0.10 

4 

514.5 

0.1651 

0.15 

6 

501.7 

0.1572 

0.15 

6 

501.7 

0.1601 

0.15 

6 

501.7 

0.1624 

R-value 


0.985 

0.982 

0.989 

0.988 

0.987 

0.986 

0.977 

0.972 

r.971 

0.970 

0.968 

0.972 

0.969 

0.965 

0.968 

0.964 

0.968 

0.965 

0.966 

0.969 

0.973 

0.973 

0.968 

0.892 

0.880 

0.946 

0.964 

0.938 

0.921 

0.963 

0.965 

0.964 

0.966 

0.949 

0.976 

0.942 

0.912 

0.949 

0.993 

0.991 

0.995 

0.995 

0.993 

0.992 

0.992 

0.992 
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-A4>ha 


R-value 


P 


Incxement  Total  Distance 


Irnm) _ (mm) 


0.15 

6 

0.15 

6 

0.15 

6 

0.10 

5 

0.10 

5 

0.10 

5 

0.10 

5 

0.10 

5 

0.10 

5 

0.15 

6 

0.15 

6 

0.15 

6 

0.15 

6 

0.15 

6 

0.15 

6 

0.15 

6 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

3 

0.10 

3 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.10 

4 

0.50 

20 

0.25 

15 

0.25 

10 

0.50 

20 

0.25 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

0.50 

20 

Wavelength 

_ (om) _ 


496.5 

0.1656 

496.5 

0.1637 

496.5 

0.1631 

488.0 

0.1775 

488.0 

0.1782 

488.0 

0.1794 

476.5 

0.1682 

476.5 

0.1701 

476.5 

0.1640 

472.4 

0.1784 

472.4 

0.1473 

472.4 

0.1612 

472.4 

0.1615 

465.8 

0.1517 

465.8 

0.1469 

465.8 

0.1556 

457.9 

0.1154 

457.9 

0.1086 

457.9 

0.0882 

454.5 

0.1776 

454.5 

0.1746 

454.5 

0.1752 

472.7 

0.0589 

476.5 

0.0629 

476.5 

0.0863 

476.5 

0.0749 

488.0 

0.0851 

488.0 

0.0758 

488.0 

0.0925 

488.0 

0.0817 

528.7 

0.0753 

528.7 

0.0770 

700.0 

0.0767 

700.0 

0.1063 

710.0 

0.1006 

710.0 

0.0962 

710.0 

0.1097 

725.0 

0.0953 

725.0 

0.0979 

725.0 

0.0977 

733.0 

0.1034 

733.0 

0.1014 

733.0 

0.1019 

750.0 

0.1201 

750.0 

0.1162 

750.0 

0.1159 

765.0 

0.1156 

765.0 

0.1140 

765.0 

0.1140 
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0.989 

0.988 

0.989 

0.995 

0.996 

0.996 

0.995 

0.996 

0.995 

0.992 

0.991 

0.993 

0.990 

0.990 

0.990 

0.992 

0.973 

0.956 

0.955 

0.992 

0.994 

0.994 

0.951 

0.976 

0.974 

0.952 

0.965 

0.967 

0.965 

0.952 

0.976 

0.989 

0.982 

0.972 

0.995 

0.998 

0.996 

0.995 

0.995 

0.994 

0.991 

0.992 

0.989 

0.968 

0.972 

0.973 

0.978 

0.980 

0.977 


Increment 
_ (mml _ 

Total  Distance 
rrem) 

Wavelength 

fnm) 

-A^)ha 

R-value 

20 

775.0 

0.1187 

0.974 

0.50 

20 

775.0 

0.1150 

0.976 

0.50 

20 

775.0 

0.1174 

0.973 

0.50 

20 

785.0 

0.1120 

0.973 

0.50 

20 

785.0 

0.1144 

0.972 

0.50 

20 

785.0 

0.1123 

0.976 

0.50 

20 

800.0 

0.1060 

0.968 

0.50 

20 

800.0 

0.1074 

0.969 

0.50 

20 

800.0 

0.1088 

0.968 

0.50 

20 

815.0 

0.1002 

0.957 

0.50 

20 

815.0 

0.1001 

0.956 

0.50 

20 

815.0 

0.1007 

0.958 

0.50 

30 

863.0 

0.1310 

0.997 

0.50 

30 

863.0 

0.1247 

0.990 

0.65 

30 

875.0 

0.1292 

0.996 

0.50 

25 

875.0 

0.1238 

0.995 

0.50 

25 

885.0 

0.1189 

0.999 

0.50 

25 

885.0 

0.1426 

0.998 

0.50 

25 

900.0 

0.1941 

0.996 

0.50 

25 

900.0 

0.1909 

0.996 

0.50 

25 

915.0 

0.1864 

0.995 

0.50 

25 

915.0 

0.1810 

0.995 

0.50 

25 

925.0 

0.1778 

0.996 

0.50 

25 

925.0 

0.1749 

0.993 

0.50 

25 

933.0 

0.1768 

0.995 

0.50 

25 

933.0 

0.1743 

0.996 

0.50 

25 

950.0 

0.1739 

0.994 

0.50 

25 

950.0 

0.1771 

0.995 

0.50 

25 

965.0 

0.1698 

0.994 

0.50 

25 

965.0 

0.1789 

0.996 

0.50 

25 

971.7 

0.1623 

0.994 

0.50 

25 

971.7 

0.1791 

0.994 

0.50 

25 

975.0 

0.1696 

0.997 

0.50 

25 

975.0 

0.1793 

0.988 

of-plane  scattering  method.  Second  3"  wafer 

Increment 

Total  Distance 

Wavelength 

-A4)ha 

R-value 

0.50 

45 

900.0 

0.992 

0.50 

45 

900.0 

0.989 

0.50 

40 

915.0 

0.992 

0.50 

40 

915.0 

0.985 

0.50 

40 

915.0 

0.992 

0.50 

40 

924.5 

0.1246 

0.994 

0.50 

40 

935.0 

0.1408 

0.988 
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Increment 

fmm) 

Total  Distance 

Wavelength 

-A^ha 

R-value 

0.50 

45 

935.0 

0.991 

0.50 

40 

945.0 

0.1176 

0.997 

0.50 

40 

945.0 

0.1017 

0.991 

0.50 

40 

945.0 

0.1061 

0.992 

0.50 

40 

885.0 

0.0727 

0.995 

0.50 

40 

885.0 

0.1099 

0.988 

0.50 

40 

875.0 

0.1083 

0.991 

0.50 

40 

875.0 

0.1289 

0.985 

0.50 

45 

850.0 

0.1832 

0.986 

0.50 

45 

850.0 

0.1826 

0.990 

0.50 

40 

850.0 

0.1950 

0.990 

0.50 

40 

850.0 

0.1836 

0.991 

0.50 

40 

850.0 

0.1512 

0.997 

0.50 

40 

850.0 

0.1280 

0.995 

0.50 

40 

850.0 

0.0448 

0.992 

0.50 

40 

875.0 

0.1571 

0.984 

0.50 

40 

875.0 

0.1844 

0.997 

0.50 

40 

875.0 

0.1113 

0.990 

0.50 

40 

875.0 

0.1606 

0.987 

0.50 

40 

875.0 

0.2049 

0.967 

0.50 

40 

900.0 

0.1815 

0.997 

0.50 

40 

900.0 

0.1407 

0.993 

0.50 

40 

900.0 

0.1688 

0.986 

0.50 

40 

900.0 

0.1623 

0.983 

0.50 

40 

925.0 

0.0993 

0.985 

0.50 

40 

925.0 

0.1846 

0.980 

0.50 

40 

925.0 

0.0897 

0.993 

0.50 

40 

950.0 

0.1860 

0.994 

0.60 

35 

950.0 

0.1521 

0.997 

0.60 

35 

950.0 

0.0905 

0.995 

0.60 

35 

950.0 

0.1031 

0.989 

0.60 

35 

975.0 

0.1052 

0.995 

0.60 

35 

975.0 

0.0863 

0.997 

0.60 

35 

975.0 

0.1147 

0.993 

0.50 

30 

975.0 

0.2155 

0.995 

-prism  method.  Second  3" 

wafer 

Increment 

Total  Distance 

Wavelength 

-Alpha 

R-value 

0.50 

20 

975.0 

0.0735 

0.989 

0.50 

17 

975.0 

0.1012 

0.998 

0.50 

17 

950.0 

0.1062 

0.998 

0.50 

17 

950.0 

0.1082 

0.999 

0.50 

17 

950.0 

0.1171 

0.999 

0.50 

17 

925.0 

0.1099 

0.997 
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Incxement 

fmm) 

Total  Distance 

Wavelength 

-A^ha 

R-value 

0.50 

17 

925.0 

0.1132 

0.997 

0.50 

15 

900.0 

0.1116 

0.994 

0.50 

15 

900.0 

0.0997 

0.990 

0.50 

15 

900.0 

0.0915 

0.998 

Two-prism  method,  Third  2”  wafer 


Incrensent 

(mm) 

Total  Distance 

Wavelength 

fnm) 

-Alpha 

R-value 

0.50 

12 

850.0 

0.1251 

0.995 

0.50 

12 

850.0 

0.1285 

0.988 

0.50 

12 

850.0 

0.1267 

0.992 

0.50 

12 

850.0 

0.1291 

0.989 

0.50 

12 

875.0 

0.1098 

0.996 

0.50 

12 

875.0 

0.1153 

0.995 

0.50 

12 

875.0 

0.1097 

0.995 

0.50 

12 

875.0 

0.1135 

0.991 

0.50 

12 

900.0 

0.0885 

0.996 

0.50 

12 

900.0 

0.1015 

0.999 
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APPENDIX  C:  Planar Waveguide.m 


DOCUMENTATION 

To  run  the  program,  open  up  a  Mathematica  shell  and  type  on  the  first  line: 

«[•*  PlanarWaveguides"  ”  ] 

This  will  open  the  package  PlanarWaveguides.m  and  evaluate  ±e  apprc^ate  initialization 
cells.  Then  type  in  the  commands  necessary  as  outlined  below.  When  done  entering  the 
specifications,  go  to  the  right-hand-side  of  the  shell  and  click  on  the  outline.  When  the 
subshell  is  outlMed,  hit  the  <entet>  key  (not  the  <retum>).  A  box  will  probably  appear 
asking  you  to  locate  the  PlanarWaveguides.m  package  so  you  must  locate  the  package  and 
then  it  will  run  to  completion. 

The  following  are  input  guidlines  for  each  subroutine: 

PropagationConstant[nl,n2,n3,n4,cladding.guide,wavelength]  finds  all  roots  for  the 
parameters  entered  and  plots  the  waveform. 

N1  Attenuation[  ( start,end,increment}  ,n2m3,n4,cladding,guide,wavelength] 
calculates  attenuation  (in  dB/cm)  for  the  range  nl=staTt,  nl<en^  nl-t-=increment  and 
plots  the  results.  The  input  must  be  entered  in  the  outlined  format. 

N2Attenuation[n  1_,  ( start,end4ncrement}  4i3,n4, cladding, guide,  wavelength] 
calculates  attenuation  (in  dB/cm)  fOT  the  range  n2=start,  n2<end,  n2+«increment  and 
plots  the  results.  The  input  must  be  entered  in  the  outlined  format. 

N3Atienuation[n  1  ,n2,  { start,end4ncrement }  ,n4,cladding,guide,  wavelength] 
calculates  attenuatitxi  (in  dB/cm)  for  the  range  n3=start,  n3<end,  n3+=increment  and 
plots  the  results.  The  input  must  be  entered  in  the  outlined  format. 

N4Attenuation[nl,n24}3,  ( start,end,increment}  .cladding, guide,  wavelength] 
calculates  attenuation  (in  dB/cm)  for  the  range  n4=star^  n4<end,  n4+=incranent  and 
plots  the  results.  The  input  must  be  entered  in  the  outlined  framat. 

CladdingAttenuation[n  1  ,n2,n3,n4,  { start,end,increment }  .guide,  wavelength] 
calculates  attenuation  (in  dB/cm)  for  the  range  cladding=start,  cladding<end, 
cladding+=increment  and  plots  the  results.  The  input  must  be  entered  in  the  outlined 
format. 

GuideAttenuation[nl,n2,n3,n4,cladding,{start,end4ncrement},wavelength] 
calculates  attenuation  (in  dB/cm)  for  the  range  gui(fe=start,  guide<end, 
guide+=inctement  and  plots  the  results.  The  input  must  be  entered  in  the  outlined 
format. 

WavelengthAttenuation[n  1  ,n2,n3,n4,cladding, guide,  { start,end,increment }] 
calculates  attenuation  (in  dB/cm)  for  the  range  wavelength=start,  wavelen^<end, 
wavelength+=nncrement  and  plots  the  results.  The  input  must  be  entered  in  the 
outlined  formaL 


68 


CutoffThickness[{stan,end4ncrement}^1^3,wavelength]  finds  the  cutoff  thickness 
for  the  parameters  entered.  The  input  must  be  entered  in  the  outlined  format 

CutoffWavelength[{startend4nciement),nl,n3,guide]  finds  the  cutoff  wavelength 
for  the  parameters  entered.  The  input  must  be  entered  in  the  outlined  format 

SOURCE  CODE 


( *  :Title:  PlanarWaveguides  ♦) 

(*  :Author:  Angela  McPherson  *) 

( *  rSummary;  There  are  four  sections: 

1.  PropagationConstant  generates  plots  of  the  waveform 
for  all  modes  found  for  the  given  parameters. 

2.  Attenuation 

a.  N1  Attenuation  generates  values  of  attenuation 
for  a  range  of  nl  values  and  plots  the  output. 

b.  N2Attenuation  (same  format  as  above) 

c.  NSAttenuation 

d.  N4 Attenuation 

e.  CladdingAttenuation 

f.  GuideAttenuation 

g.  Wavelength  Attenuation 

3.  CutoffThickness  generates  the  cutoff  values  for  the 
guiding  layer  thickness  below  which  there  will  be  no 
confined  modes. 

4.  CutoffWaveiength  (same  format  as  above). 


( *  rContext:  PlanarWaveguides'  *) 

(*  iPackage  Version:  2.1  *) 

( *  :Key  words:  waveguide,attenuation,propagation  constant, 

mode  *) 

(*  :Source: 

Borland,  William  C.  1989.  "Modeling  Optical  Waveguide  Dispersion:  . 
Four-Layer  Planar,  Channel,  and  Quantum  Well  Structures."  Ph.  D. 
Dissertation.  Univeristy  of  Cincinnati  Department  of  Electrical  and 
Computer  Engineering. 

Yariv,  Amnon.  1985.  Optical  Electronics  (Third  Edition).  New  York,  CBS 
College  Publishing. 

Zelmon,  David  E.,  William  C.  Borland,  Howard  E.  Jackson,  Joseph  T.  Boyd. 
1986.  "Near  Cutoff  Propagation  in  Low  Loss  Optical  Waveguides  Formed 
On  Si02/Si  Substrate."  in  SPIE  Vol  704  Integrated  Optical  Circuit 
Engineering  IV. 
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'.Mathematica  Version:  2.2  *) 


*) 

(* 

(*  iLimitation;  The  length  of  time  involved  depends  on  the 

number  of  iterations  asked  for  and  how  quickly  the  solutions 
converge.  On  average,  the  calculation  takes  several  minutes. 

*) 

( *  :Discussion:  Before  using  this  package,  see  the  ::usage 
information  for  details  on  the  input  format.  In  places 
where  there  is  a  range,  that  is  the  range  of  the  parameter 
you  wish  to  vary  (as  stated  in  the  procedure  name.) 

*) 


BeginPackage["PlanarWaveguides''*] 


FindRoots::usage="\nFindRoots[increment2_.threeLayerEquation[  variable].,  fo 
urLayerEquation[variable]_,k_,n2_,n3_,listCount_,flag_.list_]  is  the  core  of  the 
code  that  scans  the  range  of  possible  roots  to  find  all  confined  mode 
parameters" 

PIotThis::usage="\nPlotThis(iist_,plotLabel_,biggcstListCount_]  is  the  plotting 
routine  for  the  attenuation  plots" 

CheckInput::usage="\nCheckInput[nl_,n2_,n3_.,n4_,halfGuide_,  wavelength., k 
_,start_,end_,increment_]  checks  to  see  if  the  input  values  satisfy  appropriate 
restraints.  Start,  end  and  increment  are  options.” 

PropagationConstant::usage="\nPropagationConstant[nl_,n2_,n3_,n4_,cladding 
_,guide.,wavelength_]  finds  all  roots  for  the  parameters  entered  and  plots  the 
waveform." 

N1  Attenuation:  :usage«"\nNlAttenuation[{start_,end_, increment.  ),n2_,n3_,n4_, 
cladding.,guide.,wavelength.]  calculates  attenuation  (in  dB/cm)  for  the 

range  nl=start.  nl<end,  nl-f=increment  and  plots  the  results.  The  input  must 
be  entered  in  the  outlined  format." 

N2Attenuation:  :usage="\nN2Attenuation[nl_,{  start., end., increment.  },n3.,n4., 
cladding.,guide.,wavelength.]  calculates  attenuation  (in  dB/cm)  for  the 

range  n2=start,  n2<end,  n2-f=increment  and  plots  the  results.  The  input  must 
be  entered  in  the  outlined  format." 

N3  Attenuation:  :usage="\nN3Attenuation(nl.,n2.,{  start., end., increment.  },n4., 
cladding.,guide.,wavelength.]  calculates  attenuation  (in  dB/cm)  for  the 

range  n3=start,  n3<end,  n3-i-=increment  and  plots  the  results.  The  input  must 
be  entered  in  the  outlined  format." 
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N4  Attenuation:  ;usagc=’*\nN4Attcnuation[nl_,n2_.n3_.{stari_,end_,  increment. } , 
cladding_.gui<le_.wavelength_]  calculates  attenuation  (in  dB/cm)  for  the 
range  n4sstait,  n4<end.  n4-i>=increment  and  plots  the  results.  The  input  must 
be  entered  in  the  outlined  format." 

CladdingAttenuation::usage=”\nCladdingAttenuation[nl_,n2_,n3_,n4_,{  start., e 
nd_.increment.), guide., wavelength.]  calculates  attenuation  (in  dB/cm)  for 
the  range  cladding=start,  cladding<end,  cladding-i-=increment  and  plots  the 
results.  The  input  must  be  entered  in  the  outlined  format.” 

GuideAttenuation;:usage="\nGuideAttenuation[nl.,n2_,n3_,n4_,cladding_,{star 
t_,cnd_,incremeni_), wavelength.]  calculates  attenuation  (in  dB/cm)  for  the 
range  guide=start,  guide<end,  guide-t-=increment  and  plots  the  results.  The 
input  must  be  entered  in  the  outlined  format.” 

WavelengthAttenuation::usage=”\nWavelengthAttenuation[nl_,n2.,n3_,n4_,cl 
adding_,guide.,{start.,cnd.,increment_}]  calculates  attenuation  (in  dB/cm)  for 
the  range  wavelength=start,  '  avelength<end,  wavelength-f-=increment  and 
plots  the  results.  The  input  must  be  entered  in  the  outlined  format." 

CutoffThickness::usage=”\nCutoffThickness(  {start.,end_, increment. }  ,nl_,n3., 

wavelength _ ]  finds  the  cutoff  thickness  for  the  parameters  entered.  The 

input  must  be  entered  in  the  outlined  format.” 

CutoffWavelength:;usagc=”\nCutoffWavelength[  {start., end., increment.), nl_, 
n3.,guide.]  finds  the  cutoff  wavelength  for  the  parameters  entered.  The  input 
must  be  entered  in  the  outlined  format.” 


Off{Greater::nord]; 
Off[Less::nord]; 
Off[General:  rspell] ; 
Off[General::speIll]; 
Off[ReplaceAll::reps]; 


Begin["'private'"] 


CheckInput[nl.,n2.,n3.,n4.,halfGuide_,  wavelength. ,k_, start. :0,end_:0,incre 
ment.:0]  := 

Block[  {cutoffW  ave,cutoffGuide } , 
cutoffWave  =  .:cutoffGuide  = 

If  [(n  1<  1  )ll(n2<  1  )ll(n3<  1  )li(Re[n41<  1 ), 

Print[''One  of  the  indices  of  refraction  entered  is  less  than  1.”]; 

Break[]; 

]; 

If[(((Re[start]>Re[end])&&(Re[increment]>=0))II 

((lm[start]>lm[end])&&(lm[increment)>=0))), 

Print["The  input  range  values  were  entered  incorrectly”]; 

Break]]; 
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]; 

If[((nl>=n2)ll(n3>=n2)), 

Prinl["The  basic  requirement  for  waveguiding  (nl<n2>n3)  was  not  satisfied"]; 
Break[]; 

]; 

cutoffGuide  =  ArcTan(Sqn[(n3''2  -  nl''2)/(n2''2  -  n3''2)]]/(2  k  Sqrt[n2'^2  - 
n3''2]); 

IfKcutoffGuide  >  halfGuide), 

Print["The  guiding  thickness  entered  is  below  cutoff']; 

Brcak[]; 

]; 

cutoffWave  =  4  halfGuide  N[Pi]  Sqrt(n2'^2  -  n3''2]/ArcTan[Sqrt[n3''2- 
nlA2]/Sqrt(n2'^2-n3''2]]; 

If[(cutoffWave  <  wavelength), 

Print["The  wavelength  entered  is  below  the  cutoff  wavelength"]; 

Break[]; 

]: 


FindRoots[increment2_.variablc_,threeLayerEquation_,fourLayerEquation_,k_ 
,n2_,n3_,flag_]  := 

Block[{  rule  l,rule2, test, newList], 
ncwList  =  Table[0,{x,0,IJ,{y, 0,1]]; 

Fur[test=k  n3,test<k  n2,test+=increment2, 
rulel  =  .;rule2  =  .; 

condition  1  =  threeLayerEquation/.B->test; 
condition2  =  threeLaycrEquation/.B->(test-increment2); 
lf[(conditionl>0  &&  condition2<0)  II 
(condition  1<0  &&  condition2>0), 

nilel  =  FindRootfthreeLayerEquation  ==  0, 

{B,tcst),  Accuracy  Goal  ->  6]; 

If[NumberQ[B/.rulel], 

Check[rulc2  = 

FindRoot[fourLayerEquation  ==  0, 

{B,B/.rulcl},MaxIterations  ->  20, 

AccuracyGoal  ->  6],Continuc[]]; 

If[(((Re[B/.rulel]>k  n3)&&(Re[B/.rule2]>k  n3))&& 

((Re[B/.rulel]<k  n2)&&(Re[B/.rulc2]<k  n2))&& 
(Im[B/.rule2]<0))ll((Re[B/.rule2]>k  n3)&& 

(Re[B/.rule2]<k  n2)&&(Im[B/.rule2]<0)). 

Iftflag=l, 

ncwList  =  Union(rewList,{{2  variable,4.34  Im[B/.rule2]  10''(-4)}}], 
ncwList  =  Union[newList,{  {variable,4.34  Im[B/.rule2]  10''(-4)]}]; 

]: 

]; 

]; 

]; 

If[(Re[B/.rule2]>test),test  =  Re(B/.rule2]]; 
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If[Length[Complement[newList.{  (0.0)  }]]~0, 

Print[”Thcre  are  no  guided  modes  for  the  four  layer  guide  at  ".variable] 

]; 

newList  =  Complement[newList.{{0,0})]; 


For(loop=2,loop<=Lcngth[newListJ,loop++, 
If[Abs[Re[newLisi[[loop,2]]-newList[(loop-l,2]]]]<=.001, 
newList  =  Delete[newList.loopl; 

1; 

]; 

For[loop=  1  ,loop<=Length[newList].loop++, 

Print[newList(  [loop]  ]] ; 

]; 

Reium[newList]; 

1: 


PlotThis[list_,plotLabel_,biggestListCount_]  := 
Block[  (newList.modeFlags,i.j,list2  ] , 
newList  =  Complement[list,{  {0.0}}]; 
modeFlags  =  Tablc[0,{x.0,Length[newList]}); 
For[j=  1  ,j<=biggcstListCount,j++, 
i=j; 

While((i<=Length[newList]), 
While((ncwList[[i,l]]  =  newList([(i+l),l]]), 


]; 

Ifl(modeFlags[[i]]!=0)&&((newList[[(i-j+l),l]]= 

ncwList[[(i-j+2),l]])ll(newList[[(i-j+3),l]]== 

newList([(i-j+2),l]]))&&(modcFlags[[i-j+l]]==0), 

modeFlags([i-j+l]]=j; 

3; 

inmodeFlags[[i]]==0,modcFlags[[i]]=j]; 

i-H-; 

]; 

]; 

For[j= 1  ,j<=biggestListCount,j++. 

Iist2  =  Table[0.{x,0,l},{y.0.1}]; 
For[i=l,i<=LengUi[newList],i++, 

If[modeFlags[[i]]==j, 

list2  =  Union[list2,{{Re[newList[[i.l]]],newList[[i,2]]}}]; 

1; 

]: 

list2  =  Complement[list2,{{0,0}}l; 

Print[list2]; 

ListPlot[list2,PlotStyle  ->  PointSize(.02], 

PlotRange  ->  All.PlotLabcl  ->  plotLabel]; 

3: 
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1; 


PropagationConstant[nl_,n2_,n3_,n4_,cladding_,gui<lc_,wavelcngth_) 
Block[  {q.h.p.s, equations. equation4,k.B. increment. list}. 
halfGuide  =  .5  guide; 
k  =  2*N[Pi]/wavelength; 

q  =  Sqrt(B''2  -  k'^2  nl''2];h  =  Sqrt[k'^2  n2'^2  -  B^2]; 
p  =  Sqrt[B'^2  -  n3'^2];s  =  Sqrt(k''2  n4''2  -  B'^2]; 

equations  =  Tan{2  halfGuide  h]  -  h  (p  +  q)/(h''2  -  p  q); 
equation4  =  ((p+I  s)  ((h^2-p  q)  Tan[2  halfGuide  h]- 
h  (p+q))+Exp[-2  p  cladding]  (p-I  s)  ((h'^2+p  q)* 

Tan[2  halfGuide  h]+h  (p-q))); 
increment  =  (k  n2  -  k  n3)/1000; 
list  =  Table[0.{x.0.1}]; 

Checklnput[nl.n2.n3.n4. halfGuide.  wavelength. k]; 
lf[increment==0.increment  =  1000]; 

For[test=k  n3.test<k  n2.test+=increment. 
rulel  =  .;rule2  =  .; 
condition  1  =  equation3/.B->test; 
condition2  =  equation3/.B->(test-increment); 

I  f[  (condition  1>0  &&  condition2<0)  II 
(condition  1<0  &&  condition2>0), 
rulel  s  FindRoot[equation3  =  O.fB.test). 

AccuracyGoal  ->  6]; 

If[NumberQ[B/.rulel]. 

Check(rule2  =  FindRoot[equation4  ~  0.{B.B/.rulel}. 

Maxlterations  ->  20.  AccuracyGoal  ->  10], Continue}]]; 
IfT(((Re[B/.rulel]>k  n3)&&(Re(B/.rule2]>k  n3))&& 

((RelB/.rulel]<k  n2)&&(Re[B/.rule2]<k  n2))&& 
(Im(B/.rule2]<0))l!((Re[B/.rule2]>k  n3)&& 

(Re[B/.rule2]<k  n2)&&(Im(B/.rule2]<0)). 
list  =  Union[list.{B/.rule2}]; 

]; 

If[(Re[B/.rule2]>test).test  =  Re[B/.rule2]]; 

]; 

]; 

]; 

list  =  Complement[iist.{0.B)]; 

For[loop=2.1oop<=Length[list],loop++. 
If[Abs{Re[list[[loop]]-list([loop-l]]]]<=.001. 
list  =  Delete[list.loop]]]; 

For[loop=l.loop<=Length[iist],loop++.Print[list[[loop]]]]; 

Print["Begin  plotting  routine . '']; 

For[i=l,i<=Length[list].i++. 

B  =  list([i]]; 

q  =  Sqrt[B'^2  -  k''2  nl'^2];h  =  Sqrt[kA2  n2''2  -  B^2]; 
p  =  Sqrt[B''2  -  k''2  n3''2];s  =  Sqrt[k'^2  n4''2  -  B''2]; 

MM=  1; 

CC  =  MM*(-(E''(cladding*p)*(-Sin(-(h*halfGuide)]'^(-l)  - 
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(Cos[-(h*halfGuidc)]*(-(p/(-((h*Cos{-(h*halfGuide)]''2)/ 
Sin[-(h*halfGuidc)])  -  h*Sin[-(h*halfGuidc)]))  + 
(h*CosI-(h*halfGuidc)])/(Sin(-(h*halfGuidc)]* 
(-((h*Cos[-(h*halfGuide)]A2)/Sin(-(h*halfGuide)])  - 
h*Sin(-(h*halfGuide)l))))/Sin[-(h*halfGuide)]))  + 
(-l/(2*E'^(cladding*p))  +  (I/2*s)/(E'^(cladding*p)*p))* 
(.Sm[-(h*halfGuide)]'^(-l)  - 

(Cos[-(h*halfGuidc)l*(p/(-((h*Cos[-(h*halfGuidc)]''2)/ 
Sin[-(h*halfGuide)])  -  h*Sin[-(h*halfGuidc)])  + 
(h*Cos[-(h*halfGuide)])/(Sin[-(h*halfGuide)]* 
(-((h*Cos[-(h*halfGuide)]'^2)/Sin[-(h*halfGuide)])  - 
h*Sin[-(h*halfGuidc)]))))/Sin(-(h*halfGuidc)]  - 
E''(2*cladding*p)*(-Sint-(h*halfGuidc)]'^(- 1)  - 
(Cos(-(h*halfGuidc)]*(-(p/(-((h*Cos[-(h*halfGuidc)r2)/ 
Sin[-(h*halfGuide)])  -  h*Siii[-(h*halfGuide)]))  + 
(h*Cos[-(h*halfGuide)])/(Sin(-(h*halfGuide)]* 
(-((h*Cos[-(h*halfGuidc)]''2)/Sin[-(h*halfGuide)])  - 
h*Sin[-(h*halfauidc)]))))/Sin[-(h*halfGuide)]))); 

DD  =  (-(E^(cladding*p)*(-(p/(-((h*Cos[-(h*halfGuidc)]'^2)/ 
Sin[-(h*halfGuide)])  -  h*Sin[.(h*halfGuide)]))  + 
(h*Cos[-(h*halfGuide)])/(Sin[.(h*halfGuidc)]* 
(-((h*Cos[-(h*halfGuide)]''2)/Sin(-(h*halfGuidc)])  - 
h*Sin[-(h*halfGuidc)l))))  +  (-l/(2*E'^(cladding*p))  + 
(I/2*s)/(E''(cladding*p)*p))*(p/ 
(-((h*Cos[-(h*halfGuidc)]''2)/Sin[-(h*halfGuidc)])  - 
h*Sin[-(h*halfGuidc)])  +  (h*Cos(-(h*halfGuide)])/ 
(Sin[-(h*halfGuidc)]*(-((h*Cos[-(h*halfGuide)]''2)/ 
Sin[-(h*halfGuide)])  •  h*Sin[-(h*halfGuide)]))  - 
E''(2*cladding*p)*(-(p/(-((h*Cos[-(h*halfGuide)]''2)/ 
Sin[-(h*halfGuide)])  -  h*Sin[-(h*halfGuidc)]))  + 
(h*Cos[-(h*halfGuide)])/(Sin(-(h*halfGuide)]* 
(-((h*Cos[-(h*halfGuide)]''2)/Sin[-(h*halfGuide)])  - 
h*Sin[.(h*halfGuide)]))))); 


FF  =  MM*(E''(cladding*p)  +  E''(2*cladding*p)* 
(-l/(2*E''(cladding*p))  +  (I/2*s)/(E'^(cladding*p)*p))); 

GG  =  MM*(l/(2*E''(cladding*p))  + 
(-I/2*s)/(E''(cladding*p)*p)); 

BB  =  CC  Sin[h  halfGuide]  +  DD  Cos[h  halfGuide]; 


plotl  =  Plot[Re[BB  Exp[-q  (x.  -  halfGuide)]], 
{x,halfGuide,2  halfGuide J.DispIayFunction  ->  Identity]; 

plot2  =  Plot[Re[CC  Sin[h  x]  +  DD  Cos[h  x]], 

{x,  -halfGuide,  halfGuide), DisplayFunction  ->  Identity]; 

ploG  =  PIot[Re[FF  Exp[p  (x  +  halfGuide)]  + 

GG  Exp[-p  (x  +  halfGuide)]], {x,  -(cladding  +  halfGuide), 
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-halfGuide). Display  Function  ->  Identity]; 


plot4  =  Plot[Re[MM  Exp[I  s  (r  +  cladding  +  halfGuide)]]. 

(x.  -1.5  (cladding  +  halfGuide).  -(cladding  +  halfGuide)}. 

DisplayFunction  ->  Identity]; 

Show[plot4.plot3,plot2,plotl.  PlotRange  ->  All. 

PlotLabel  ->  "Wave  Function". 

DisplayFunction  ->  SDisplayFunction. 

Ticks  ->  {{halfGuide.-haifGuide,-(cladding  +  halfGuide)}. 

Automatic}]; 

]; 

]  /; 

NumberQ[N[nl]]&&NumberQ[N[n2]]&&NumberQ[N(n3]]&&NumbcrQ[N[n4]]&&Nu 

mberQ[N[cladding]]&&NumberQ[N[guide]]&&Number(^[N[wavelength]] 


N1  Attenuation}  {start_,cnd_.increment_}  .n2_,n3_.n4_.cladding_,guide_.wavele 
ngthj  ;= 

Biock[{list.increment2.variable,iistCount.biggestListCount, 
plotLabel. flag. halfGuide. q.h.p.s.k.fourLayerEquation, 
threeLayerEquation.argument}. 
biggestListCount  =  0; 
flag  =  0; 

list  =  Table[0,{x.0.1}.{y.0.1}]; 
halfGuide  =  .5  guide; 

q(argument_]  :=  Sqrt(B''2  -  k''2  (argument)''2]; 
h  =  Sqrt[k'^2  n2''2  -  B''2]; 
p  =  Sqrt[B''2  -  k''2  n3''2]; 
s  =  Sqrt[k^2  n4'^2  -  B''2]; 
k  =  2  N[Pi]/wavelength; 

fourLayerEquation[argument_]  :=  ((p  +  I  s)  ((h^2  -  p  q[argument])  Tan[2 
halfGuide  h]  -  h  (p  +  q[arguraent]))  +  Exp[-2  p  cladding]  (p  -  I  s)  ((h''2  +  p 
q[argument])  Tan[2  halfGuide  h]  +  h  (p  -  q[argument]))); 

threeLayerEquation[argument_]  ;=  Tan[2  halfGuide  h]  -  h  (p 
q[argument])/(h^2  -  p  q[argument]); 

Checklnput[start. n2.n3,n4. halfGuide, wavelength.k. start. end, increment]; 
Checklnput[end,n2, n3,n4, halfGuide,  wavelength.k, start, end, increment]; 
lf[increment==0,increment  =  1000]; 
plotLabel  =  "Attenuation  (dB/cm)  vs  nl”; 
increment2  =  (k  n2  -  k  n3)/10; 
For[variable=start,variable<=end.variable-)-=increment, 
list  = 

Union  [1  is  t.(list2=:FindRoots[  increment  2,  vari  able,  threeLayerEquation[vari  able  ] , 
fourLayerEquation}  variable], k,n2,n3. flag])]; 

If[Length[Complement[list2.(  (0,0}  ]]]>biggestListCount, 
biggestListCount=Length[Complement(list2,{  {0,0} }]]; 

]; 

]; 

PlotThis[list,plotLabel, biggestListCount]; 


]  /: 

NuinberQ[N[start]]&&NuinberQ[^[end]]&&NumberQ[N[iacrement]]&&NumberQ[ 
N[n2]]&&NuinberQ[N[n3]]&&NumberQ[N[n4]]&&NumberQ[N[cladding]]&&Numb 
crQ[N[guidc]]&&NumberQ[N(  wavelength]] 


N2Aitenuation[nl_,{start_,end_,increment_],n3_.n4_,cIadding_.guide_,waveIe 
ngthj  := 

Block[{list.increment2,variable.listCount.biggestListCount, 
plotLabel.flag,halfGuide,q,h,p,s,k,fourLayerEquatioa.threeLayerEquation) , 
biggestListCount  =  0; 
flag  =  0; 

list  =  Table[0,{x.0,l}.{y.0,l}]; 
halfGuide  =  .5  guide: 
q  =  Sqn[B''2  -  k'^2  nl''2]; 

h[argument_]  :=  Sqrt(k''2  (argument)''2  -  B'‘2]; 
p  =  Sqn[B'^2  -  k'^2  n3^2]; 
s  =  Sqrtilc'^2  n4''2  -  B''2]: 
k  =  2  N[Pi]/wavelength; 

fourLayerEquation[argument_]  :=  ((p  +  I  s)  (((h[argument])'^2  -  p  q)  Tan[2 
halfGuide  h(argument]]  -  h[argument]  (p  +  q))  +  Exp[-2  p  cladding]  (p  -  I  s) 
(((h[argument])^2  p  q)  Tan[2  halfGuide  h[argument]]  +  h[argument]  (p  • 

q))); 

threeLayerEquation[argument_]  :=  Tan[2  halfGuide  h[arguinent]]  • 
h(argument]  (p  +  q)/((h[argument])'^2  -  p  q); 

Checklnput[nl, start, n3,n4,haifGuide,wavelength,k, start, end, increment]; 
Checklnput[n  1, end, n3,n4, halfGuide.  wavelength,k,start,end, increment]; 
If[increment=:sO,increment  =  1000]; 
plotLabel  =  "Attenuation  (dB/cm)  vs  n2”; 
For[variable=start,variable<=end.variable+=increment, 
increment2  =  (k  variable  •  k  n3)/10: 
list  = 

Union  [list,(list2=FindRoots[increment2,variabie.threeLayerEquation  [variable], 
fourLayerEquation[variable],k,variable,n3,flag])]; 

If[Length[Complement[list2.{  {0,0}  }]]>biggestListCount, 
biggestListCount=Length[Complement[list2,{  [0,0]  ]]]; 

]; 

]; 

PlotThis[list,plotLabel,biggestListCount]; 

]  /; 

NumberQ[N[nl]]&&NumberQ[N[stari]]&&NumberQ[N[end]]&&NumberQ[N[incre 
ment]]&&NumberQ[N[n3]]&&NumberQ[N[n4]]&&NumberQ[N[cladding]]&&Numb 
crQ[N[guide]]&&NumberQ[N[  wavelength]] 


N2ComplexAttenuation(nl_,{start_,cnd_,increment_}  ,n3_,n4_,cladding_,guide_, 
wavelength_]  := 

Block[  (list, increment2,variable.listCount, biggestListCount, 
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plotLabel,flag,halfGuide,q,h.p,s,k,fourLayerEquation,threeLayerEquation,iinag 
inary  Var.realVar) . 
biggestListCount  =  0; 
flag  =  0; 

list  =  Tablcl0.{x.0.1).{y.0.1}l: 
halfGuide  =  .5  guide; 
q  =  Sqrt[B''2  -  k''2 

h[argunient_]  :=  Sqrt[k''2  (argument)^2  -  B''2]; 
p  =  SqrttB'^2  -  k-'2  n3'^2]; 
s  =  Sqrt[k''2  n4>^2  -  B'^2]; 
k  =  2  N[Pi]/wavelcngth; 

fourLayerEquation[argument_]  :=  ((p  +  I  s)  (((h[argumcnt])''2  -  p  q)  Tan[2 
halfGuide  h[argument]]  -  h[argument]  (p  +  q))  -t-  Exp[-2  p  cladding]  (p  •  I  s) 
(((h[argument])^2  -t-  p  q)  Tan[2  halfGuide  h[argument]]  +  h[argument]  (p  - 
q))): 

threeLayerEquation [argument.]  ;=  Tan[2  halfGuide  h [argument]]  • 
h[argument]  (p  +  q)/((h[argument])'^2  -  p  q); 

Checklnput[n  1. start, n3.n4. halfGuide. wavelength, k, start, end, increment]; 

Checklnput[n  1. end. n3.n4, halfGuide.  wavelength.k.start.end. increment]; 
lf[increment==0,increment  =  1000]; 
plotl.abel  =  "Attenuation  (dB/cm)  vs  Re[n2]"; 


increment2  =  (k  n2  -  k  n3)/10; 


If[Re[increment]ssO, 
inciement2  =  (k  n2  -  k  n3)/10; 

For[imaginaryVar=Abs[Im[start]],imaginaryVar<=Abs[Im[end]],imaginaryVar-t- 

=Abs[Im[increment]]. 

variable  =  Re[start]  -  imaginaryVar  I; 
list  = 

Union  [list. (list2=FindRoots[increment2.  variable, threeLayerEquation[  variable], 
fourLayerEquation[variabie],k.n2.n3,flag])]; 

If[Length[Complement[lis^.{  {0,0}  }]]>biggestListCount, 
biggestListCount=Length[Complement[list2,[  [0,0]  ]]]; 

]; 

]; 

PlotThis[list,plotLabel,biggestListCount]; 

]; 

lf[lm[increment]==0, 

For[realVar=Abs[Re[start]], real  Var<=Abs[Re[cnd]], real  Var+=Abs[Re[increment]] 
♦ 

variable  =  realVar  -  Abs[Im[start]]  I; 
increment2  =  (k  rcalVar  -  k  n3)/10; 
list  s 

Union  [list, (list2=FindRoots[increment2. variable, threeLayerEquation[variable], 
fourLayerEquation[  variable], k,n2.n3, flag])]; 

If[Length[Complement[list2.(  (0.0)  }]]>biggestListCount. 
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biggestListCount-Length[Compiement[list2,{  {0.0} )]]; 

1: 

]; 

PlotThis[list,plotLabel.biggestListCount]; 

]; 

If[(Im[increment]!=0)&&(Re[incremem]!=0). 

lf[(Re(start]-Re[end])/Re[increiQent]>=(Im[start]-Im[end])/Im[increment]. 

For[realVap=Re[start],realVar<=Re[end],realVar-t-=Re[increinem], 

For[imaginaryVar=Abs[Im{start]].imaginaryVar<sAbs[Im(end]],imaginaryVar-t- 

=Abs[Im[increment]], 

variable  =  realVar  -  imaginaiyVar  I; 
incremeat2  =  (k  realVar  •  k  nSVlO; 
list  = 

Union[list,(list2=FindRoots[increment2. variable. threeLayerEquation[variable]. 
fourLayerEquation[variable].k.n2.n3.flag])]; 

If[Length[Complement[list2.(  {0.0}  }]]>biggestListCount. 
biggestListCount=Length{Coniplement{list2.{  {0.0} }]}; 

]; 

]; 

1; 

]; 

If{(Re{start}-Re{end})/Re{increinent}<(Im{start}-Im{end})/Im{increment], 

For[imagiaaryVar=Abs[lm[starcJJ,imagiaaryVar<sAbs[Im[end]],imagiaaryVir+ 

=Abs{Im{increment}}. 

For{realVar=Re(start].realVar<sRe{end].realVar-t'=Re{incremeDt]. 
variable  =  realVar  •  imaginaiyVar  I; 
increment2  =  (k  realVar  •  k  n3)/10; 
list  = 

Union  {list. (iist2=FindRoots{increment2.variable.threeLayerEquation{variable}. 
fourLayerEquation{  variable  j.k.n2.n3. flag})}; 

If{Length{Complement{list2,{  {0.0} }  ]]>biggestListCount. 
biggestListCoum=Lengtb{Complement{list2.{  {0,0} }]]; 

}; 

]; 

]; 

]; 

PlotThis{list,plotLabel.biggestListCount}; 

]; 

]  /; 

NumberQ{N{nl}}&&NumberQ{N{n2}]&&NumberQ{N{n3}]&&NumberQ{N{start}]&& 
NumberQ{N{end}}&&NumberQ{N{increment}}&&NumberQ{N{cladding}i&&Numb 
erQ{N{guide}]&&NumberQ{N{  wavelength]] 


N3Attenuation{nl_,n2_,{start_,end_,increment_},n4_,cladding_,guide_,wavele 
ngth_}  := 

Block{{list,increment2,variable.listCount,biggestListCount, 
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plotLabel.flag.halfGuide.q.h.p.s.k.fourLayerEquation.threeLayerEquation) , 
biggestListCount  =  0; 
flag  =  0; 

list  =  Table[0,{x,0.1}.{y,0.11]; 
halfGuide  =  .5  guide; 

q  =  SqrttB'^2  -  k''2  nl'^2];h  =  Sqit[k'^2  02^2  -  B^2]; 
p[argument_]  :=  Sqrt(B''2  -  k'^2  (argument)^2]; 
s  =  Sqrt[k'^2  n4'^2  -  B''2]; 
k  s  2  N[Pi]/wavelength; 

fourLayerEquation(argument_}  ;=  ((p[argument]  -f  I  s)  ((h^2  -  p[argument]  q) 
Tan[2  halfGuide  h]  •  h  (p[argumentj  +  q))  -f  Exp[-2  p[arguinent]  cladding] 
(p[argument]  -  Is)  ((h^2  +  p[argument]  q)  Tan[2  halfGuide  h]  +  h  (p[argument] 
•  q))); 

threeLayerEquation(arguinent_]  ;=  Tan[2  halfGuide  h]  -  h  (p[argument]  + 
q)/(b^2  -  p[argument]  q); 

Checklnput(n  l.n2, start. n4. halfGuide. wavelength.k.start.end, increment]; 
Checklnput[nl,n2.end,n4, halfGuide,  wavelength, k. start, end, increment]; 
lf[increment==0.increment  =  1000]; 
plotLabel  =  "Attenuation  (dB/cm)  vs  n3"; 
For[variable=start,variable<=end.variable+=increment. 
increment2  =  (k  n2  •  k  variable)/10; 
list  = 

Union  [list, (list2=FindRoots[increment2.variable,threeLayerEquation  [variable], 
fourLayerEquation[  variable], k,n2. variable. flag])]; 

If[Length[CompIement[list2,{  (O.O)  }]]>biggestListCount, 
biggestListCount=Length[Complement[list2,(  (0,0) }]]; 

]; 

]; 

PlotThis[list,plotLabel.biggestListCount]; 

1  /; 

NumberQ[N[nl]]&&NumberQ[N[n2]]&&NumberQ[N[start]]&&NumberQ[N[end]]& 

&NumberQ[N[increment]]&&NumberQ[N[n4]]&&NumberQ[N[cladding]]&&Numb 

erQ[N[guide]]&&NumberQ[N[wavelength]] 


N4Attenuation[nl_,n2_,n3_,{start_,end_,increment_],cladding_.guidc_,wavele 
ngthj  := 

Block[(list,increment2,variable.listCount. biggestListCount, 
plotLabel. flag, halfGuide, q,h.p,s,k,fourLayerEquation,threeLayerEquation,imag 
inary  Var.realVar] , 
biggestListCount  =  0; 
flag  =  0; 

list  =  Table[0.(x.0.1}.(y.0.1}]; 
halfGuide  =  .5  guide; 

q  =  Sqrt[B''2  -  k''2  nl''2];h  =  Sqrt[k''2  *\2''2  -  B^2]; 
p  =  Sqrt[B''2  -  k'^2  n3''2]; 
s[argument_]  :=  Sqrt[k''2  (argument)^  2  -  B''2]; 
k  =  2  N[Pi]/wavelength; 

fourLayerEquation[argument_]  :=  ((p  +  1  s[argument])  ((h^2  -  p  q)  Tan[2 
halfGuide  h]  -  h  (p  -t-  q))  +  Exp[-2  p  cladding]  (p  -  I  s[argument])  ((h^2  +  p  q) 
Tan[2  halfGuide  h]  h  (p  -  q))); 
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threeLayeiEquation[argumentJ  :=  Tan(2  halfGuide  h]  -  b  (p  q)/(h^2  -  p  q); 
Checklnput[n  l.n2.n3.start,halfGuide. wavelength, k, start, end. increment]; 
CheckInput[nl.n2,n3.end,halfGuide.wavelength,k, start, end. increment]; 
IflincrementssO.increment  s  1000]: 
plotLabel  »  "Attenuation  (dB/cm)  vs  Re[n4]”; 
increment2  =  (k  n2  -  k  n3)/10; 

If[Re[increment]==0. 

For[imaginary  Var=Abs(Im[start]].imaginary  Var<=Abs[Im[end]], imaginary  Var+ 
=Abs(Im[increment]], 

variable  =  Re[start]  -  imaginaiyVar  I; 
list  = 

Union[list,(list2=sFindRoots[increment2,variabie.threeLayerEquation[variable], 

fourLayerEquation[variablej.k.n2.n3.fIag])]; 

If[Length[Complement[list2.{  {0,0)  }]]>biggestListCount. 
biggestListCount=Lengtb[Complement{list2,{  (0,0) )]]; 

]; 

]; 

PlotThis[iist.plotLabel,biggestListCount]; 

]; 

lf[lm[increment]==0. 

For(  real  Var«Abs[Re[start]], real  Var<=Abs{Re[cnd]], real  Var+=Abs[Re[incremcnt]] 

variable  =  realVar  -  Abs[Im[start]]  I; 
list  = 

Union[list.(list2=FindRoots[increment2,variable,tbreeLayerEquation[variable]. 

fourLayerEquation{variable],k.n2,n3,flag])]; 

If{Length[Complement[list2.(  {0.0}  }]]>biggestListCount, 
biggestListCount=Length[Complement{list2,{  {0.0} }]]; 

]; 

]; 

PlotThis[list,plotLabei,biggestListCount]; 

]; 

If[(Im[increment]!=0)&&(Re[increment]!=0), 

If[(Re[start]-Re{end])/Re[increment]>=:(Im[start]-Im[end])/Im[increment], 

For[realVar=Re[start],realVar<=5Re[end],realVar+=Re[increment], 

For[imaginaryVar=Abs[lm[start]],imaginaryVar<=Abs[Im[end]].imaginaryVar-t- 

=:Abs[Im[increment]], 

variable  =  realVar  -  imaginaryVar  I; 
list  = 

Union  [list, (list2=FindRoots[increment2, variable. threeLayerEquation[variable], 
fourLayerEquation[  variable], k,n2,n3, flag])]; 

If[Lengtb[Complement[list2,  { (0,0  ]  ]  ]]>biggestListCount, 
biggestListCount=Length[Complement[list2,{  {0,0} }]]; 

]: 

]; 

]; 

]: 

If{(Re{start]-Re{end])/Re{increment]<(Im{start]-Im{end])/Im{increment], 
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For[imaginaryVar=Abs[Im[start]].imaginaryVar<=:Abs[Im[end]].imaginaryVar-«- 

=Abs[Im[increment]]. 

For[realVar=Re[start].realVar<=:Re[end].realVar+=Re[increment], 
variable  =  realVar  -  imaginaiyVar  I; 
list  = 

Union  [list. (list2=:FindRoots[increment2.  variable, threeLayerEquation  [variable], 
fourLayerEquation[variablej.k.n2.n3.flag])]; 

If[Length[Complement[list2.{  {0,0}  }]]>biggestListCount. 
biggestListCount=Length[Coinplement[list2.{  [0,0) }  ]]; 

1; 

1: 

]; 

]; 

PlotThis[list,plotLabel.biggestListCount]; 

]; 

]  /; 

NumberQ[N[nl]]&&NumberQ[N[n2]]&&NumberQ[N[n3]]&&NumberQ[N[start]]&& 

NumberQ[N[end]]&&NuinberQ[N[increment]]&&NuinberQ[N[cladding]]&&Numb 

erQ[N[guide]]&&NumberQ[N[wavelength]] 


CladdingAttenuation[nl_,n2_,n3_,n4_,{start_,end_,increment_}  ,guide_,wavele 
ngthj  ;= 

Block[[list.increment2,variable,listCount,biggestListCount. 
plotLabel,flag,halfGuide.q,h,p,s,k,fourLayerEquation,threeLayerEquation), 
biggestListCount  =  0; 
flag  =  0; 

list  =  Table[0,{x.0,l},[y,0,l}]; 
halfGuide  =  .5  guide; 

q  =  Sqrt[B''2  -  k''2  nl''2];h  =  Sqrt[k''2  n2'^2  -  B''2]; 
p  =  Sqrt[B''2  -  k''2  n3''2];s  =  Sqrt[k''2  n4''2  -  B''2]; 
k  =  2*N[Pi]/wavelengA; 

fourLayerEquation[argumentJ  :=  ((p  +  I  s)  ((h'^2  -  p  q)  Tan[2  halfGuide  h]  -  h 
(P  +  q))  +  Exp[-2  p  argument]  (p  -  I  s)  ((h'^2  +  p  q)  Tan[2  halfGuide  h]  +  h  (p  - 

q))): 

threeLayerEquation[argumentJ  :=  Tan[2  halfGuide  h]  -  h  (p  +  q)/(h''2  -  p  q); 
Checklnput[n  I,n2,n3,n4, halfGuide, wavelength, k, start, end, increment]; 
lf[increment==0,increment  =  1000]; 

plotLabel  =  "  Attenuation  (dB/cm)\n  vsVnCladding  Thickness 

(microns)"; 

increment2  =  (k  n2  -  k  n3)/10; 
For[variable=start,variable<=end,variable-»-=increment, 
list  = 

Union  [list, (list2=FindRoois[incrcmcnt2, variable, thrccLayerEquation[variablc]. 
fourLayerEquation[variable],k,n2,n3,flag])]; 

If[Length[Complement[list2,{  {0,0}  }]]>biggestListCount, 
biggestListCount=Length[Complement[list2,{  {0.0] }]]; 

]; 

]; 

PlotThis[list,plotLabel,biggestListCount]; 
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]  /: 

NumberQ[N[nl]]&&NumberQ[N[n2]]&&NumberQ[N[n3]]&&NuniberQ[N[n4]]&&Nu 

mberQ[N[start]]&&NumberQ[N[end]]&&NumberQ[N[increment]]&&NuinberQ[N[g 

uide]]&&NumberQ[N[waveleagth]l 


GuidcAttenuation[nl_,ii2_,n3_,n4_,cladding_,(start_,end_,incrcmem_}  ,wav«\c 
ngthj  :* 

Block[{list.increment2.variabie,listCount.biggestListCount. 
plotLabel.flag.halfGuide.q.h.p.s.k.fourLayerEquation.threeLayerEquation), 
biggestListCount  =  0; 
flag  =  1; 

list  =  Tablc[0.{x.0.1).{y.0.1}]; 

newStan  =  .5  start;newEnd  =  .5  end;newlncrement  =  .5  increment; 
halfGuide  =  .5  guide; 

q  =  Sqrt(B'^2  -  nl^2];h  =  Sqrt[k'^2  n2'^2  -  B'^2]; 

p  =  SqrtlB^a  -  k''2  n3''2];s  =  Sqit(k''2  n4A2  -  B'^2]; 
k  =  2*N(Pi]/wavelength; 

fourLayer^uation[argument_]  :=  ((p  +  I  s)  ((h^2  -  p  q)  Tan[2  argument  h]  -  h 
(p  +  <l))  +  Exp(*2  p  cladding]  (p  -  I  s)  ((h'^2  +  p  q)  Tan[2  argument  h]  +  h  (p  - 

q))): 

threeLayerEquation[argument_]  ;=  Tan[2  argument  h]  -  h  (p  +  q)/(h'^2  -  p  q); 
CheckInput[nl.n2,n3,n4,newStart.wavelength,k,start,end, increment]; 
CheckInput[nl,n2,n3,n4,newEnd,  wavelength, k,start, end. increment]; 
If[incrementssO,incremcnt  =  1000]; 

plotLabel  s  "  Attenuation  (dB/cm)\n  vsNnGuiding  Thickness  (microns)"; 

increment2  =  (k  n2  -  k  n3)/10; 

For[variable=newStart,variable<snewEnd,variable-t-snewIncrement, 
list  = 

Union  (list. (list2=FindRoots[increment2,variable,threeLayerEquation[variable], 
fourLayerEquation(variablej,k,n2.n3,flag])]; 

If(Length[Complement[list2.{  (0,0}  )]]>biggestListCount, 
biggestListCount=Length(Complement(list2,{  (0.0] )]]; 

]; 

1; 

PlotThisdist.plotLabel.biggestListCount]; 

]  /: 

NumberQ(N(nl]]&&NumberQ[N(n2]]&&NumbcrQ[N[n3]]&&NumbcrQ(N(n4]]&&Nu 
mberQ[N(cladding]]&&NumberQ[N[start]]&&NumberQ[N[end]]&&NumberQ[N[inc 
rcment]]&&NumberQ(N[  wavelength]] 


WavelengihAttenuation[nl_,n2_,n3_,n4_.cladding_,guide_,  {start_,end_,increm 
ent_]]  := 

Block[(  list, increment2, variable, listCount, biggestListCount, 
plotLabel, flag, halfGuide, q.h.p.s.k.fourLayerEquation.threeLayerEquation), 
biggestListCount  =  0; 
flag  =  0; 

list  =  Table(0.(x.0,l),{y,0,l}]; 
halfGuide  =  .5  guide; 
k(argument_]  :=  2  N(Pi]/argument; 
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qlargumentj  :=  Sqrt[B''2  -  (k[argument])^2  nK2]; 
h[argumentj  :=  Sqrt[(k[argumentJ)^2  n2'^2  -  B'^2]; 
p[argumeatj  :=  Sqrt[B^2  -  (k[argxunent])^2  n3''2]; 
s[argumem_]  :=  Sqrt[(k[argumcnt])^2  n4'^2  -  B''2]; 
fourLayerEquation[argument_]  :=  ((p[argument]  +  I  sfargument]) 
(((h[arguinent])^2  -  p[argument]  q[argument])  Tan[2  halfGuide  h[argumem]]  - 
h[argumem]  (p[arguinent]  +  q[argument]))  +  Exp[-2  p[argument]  cladding] 
(p[argument]  -  I  s[argumeat])  (((h[argument])^2  +  p[argument]  qtargument]) 
Tan[2  halfGuide  h[argunient]]  +  h(argument]  (p[argument]  -  q[arguinent}))); 

threeLayerEquation[argument_]  :=  Tan[2  halfGuide  h[argument]]  - 
h[argument]  (p[argument]  +  q[argument])/((h[argument])^2  -  p[argument] 
q[argument]); 

CheckInput[nl,Q2.n3,n4. halfGuide, start.k  [start], start, end.increment]; 
CheckInput[nl,n2,n3,n4,halfGuide,end,k[end]. start, end, increment]; 
lf[increment~0,increment  =  1.000]; 

plotLabel  =  "Attenuation  (dB/cm)\n  vsVnWavelength  (microns)"; 

For[variable=start,variable<=end,variable-f=:increment, 
incremcna  =  (k[variablc]  n2  -  k[variablc]  n3)/10; 
list  s 

Union[list,(list2=FindRoots[increment2,variable,threeLayerEquation[variable], 
fourLayerEquation[variable],k[variable].n2.n3,flag])]; 
If[Length[Complement[list2,(  {0,0}  }]]>biggestListCount, 
biggestListCountsLength[Compiement[list2,{  (0,0) }]]; 

]t 

1; 

PlotThis[list,plotLabel,biggestListCount]; 

]  /; 

NumberQ[N(nl]]&&NumbcrQ[N[n2]]&&NumberQ[N[n3]]&&NumberQ[N[n4]]&&Nu 

mberQ[N[cladding]]&&NumbcrQ(N[guide]]&&NumbcrQ[N[siart]]&&NumberQ[N[c 

nd]]&&NumberQ[N[increment]] 


CutoffThickncss[{start_,cnd_,incrcmcnt_},nl_,n3_,wavclcngth_]  := 

Block}  {k,list,n2,cutoffGuide } , 
k  =  2*N[Pi]/wavclcngth; 
list  =  Tablc[0  x,{x,0,l},{y,0.1}]; 

If[nl==n3, 

Print["Thc  parameters  entered  define  a  symmetric  guide  therefore  there  is 
no  cutoff"]; 

Break}]; 

]; 

lf}((start>end)&&(increment>0)), 

Print}"Input  values  were  entered  incorrectly"]; 

Break}]; 

1; 

If}((nl>=sian)&&(siari<=n3))ll((nl>=end)&&(cnd<=n3)), 

Print}  "The  parameters  entered  disobey  the  requirement  for  waveguiding 
nl<n2>n3"]; 

Break}]; 

]; 

For}n2=start,n2<=end,n2+=increment, 
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cutoffGuidc=ArcTan[Sqrt[(n3''2  -  nl'^2)/(n2>'2  -  n3''2)]]/(k  Sqrt[n2^2  -  n3'^2]); 
list=Union[list.  { {n2,cutoffGuide} }  ]; 

Print(n2,"  ’’.cutoffGuidc]; 

1; 

list  =  Complement[list,{{0.0)}]; 

ListPlot[list,PlotStyle  ->  PointSize[.02]. 

PlotLabel  ->  "Cutoff  Thickness  (microns)  vs  n2”]; 

]  /; 

NumberQ[N[start]]&&NumberQ[N[end]l&&NumberQ[N[increment]]&&NumberQ[ 

N[nl]]&&NumbcrQ[N[n3]l&&NumbcrQ[N[wavelength]] 


CutoffWavelength[  { start_.end_,incrcmcnt_}  .n  l_,n3_,guidc_]  := 

Block[  {n2,list,cutoffWave } . 
list  =  Table[0  x,{x.0,l}.{y.0.ll]; 
halfGuide  =  .5  guide: 

If[nl~n3.Print["The  parameters  entered  define  a  symmetric  guide  therefore 
there  is  no  cutofr];Break[]]; 

If[((stan>end)&&(increment>0)),Print["lnput  values  were  entered 
incorrectly  "];Break[]]: 

If[(nl>=start)  II  (start<=n3)  II  (nl>=end)  II  (end<=n3). 

Print["The  parameters  entered  disobey  the  requirement  for  waveguiding 
(nl<n2>n3)"]: 

Break[]: 

]; 

For(n2=start,n2<=end,n2+=increment. 
cutoffWave=4  halfGuide  N[Pi]  Sqrt[n2''2  -  n3'^2]/ 
ArcTantSqn[n3''2-nl''2]/Sqrt[n2''2-n3''2]]; 
list=Union[iist.(  (n2.cutoffWave} }]; 

Print[n2,"  ".cutoffWave]; 

]; 

ListPlot[list.PIotStyle  ->  PointSize[.02], 

PlotLabel  ->  "Cutoff  Wavelength  (microns)  vs  n2"]; 

]  /: 

NumbcrQ[N[start]]&&NumberQ(Nlend]l&&NumberQ[N[incremcnt]]&&NumberQ[ 

N[nl]]&&NumberQ[N[n3]]&&NumberQ[N[guide]] 


End[  (*  "'private'"  *)  ] 

Protect  [PropagationConstant.Nl  Attenuation, N2Attenuation, 

N3  Attenuation, N4Attenuation.CladdingAttenuation, 

Guide  Attenuation,  WavelengthAttenuation,FindRoots,CheckInput, 
PlotThis,CutoffThickness,CutoffWavelength] 

EndPackage[] 


END 
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